Electron Phonon Interaction and Strong 
Correlations in High-Temperature 
Superconductors: One can not avoid unavoidable 



Abstract. The important role of the electron-phonon interaction (EPI) in explaining the properties 
of the normal state and pairing mechanism in high-T^ superconductors (HTSC) is discussed. A 
number of experimental results are analyzed such as: dynamical conductivity, Raman scattering, 
neutron scattering, ARPES, tunnelling measurements, isotope effect and etc. They give convincing 
evidence that the EPI is strong and dominantly contributes to pairing in HTSC oxides. It is argued 
that strong electronic correlations in conjunction with the pronounced (in relatively weakly screened 
materials) EPI are unavoidable ingredients for the microscopic theory of pairing in HTSC oxides. 
1 present the well defined and controllable theory of strong correlations and the EPI. It is shown 
that strong correlations give rise to the pronounced /orwarii scattering peak in the EPI - the ESP 
theory. The ESP theory explains in a consistent way several (crucial) puzzles such as much smaller 
transport coupling constant than the pairing one (Ar^ ^ A), which are present if one interprets the 
results in HTSC oxides by the old Migdal-Eliashberg theory for the EPI. The ARPES shift puzzle 
where the nodal kink at 70 meV is unshifted in the superconducting state, while the anti-nodal one at 
40 meV is shifted can be explained at present only by the ESP theory. A number of other interesting 
predictions of the ESP theory are also discussed. 
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1. INTRODUCTION 
1.1. Importance of strong electronic correlations and EPI 



Seventeen years after the discovery of the high-Tc superconductors (HTSC) [1] there 
is still no consensus about the pairing mechanism in these materials. At present two 



possible theories are in the focus, the first one based on the electron-phonon interaction 
(EPI) and the second one based on spin fluctuation interaction (SFI). In the meantime it 
was well established that metallic compounds of HTSC oxides are obtained from insu- 
lating parent compounds by doping with small number of carriers - usually called holes. 
It turns out that the parent insulating state is far from being conventional band insulator 
where usually an even number of electrons (holes) per lattice site fill Bloch bands com- 
pletely. By counting the electron number one comes (naively) to the conclusion that the 
parent compounds of copper oxides (for instance LaiCuO^ and YBa2Cu^0()) should be 
metallic, because in the unit cell there is odd (nine) number of J-electrons per Cm^+ ion. 
The way out from this controversy is in the presence of strong electronic correlations. 
They are due to the localized <i-orbital on the Cm^"*" ion giving rise to the strong Coulomb 
repulsion U of two ?>d^2_y^ electrons (or holes) at a given lattice site with opposite spins. 
This repulsion keeps electrons apart making them to be localized on the lattice, but with 
localized spins (5 = 1/2). This type of insulating state is called the Mott-Hubbard insu- 
lator. Speaking in language of electronic bands, for large on-site repulsion f/ ^ W and 
for one electron per lattice site the original conduction band (with the width W) is split 
into the lower Hubbard band with localized spins and the empty upper band separated 
by U from the lower one - see more in [2] and Section 4. 

The relevance of strong correlations is well documented experimentally: (i) The 
electron-energy-loss spectroscopy [24] shows a transfer of intensity (which is a mea- 
sure of the number of states) from higher to lower energies by doping. Such a property 
is characteristic for the class of Hubbard models where the number of states in the up- 
per Hubbard band decreases by increasing the hole doping. For comparison, in typical 
semiconductors the number of states in the valence band is determined by the number 
of atoms, i.e. it is fixed and doping independent, (ii) The self-consistent band-structure 
calculations and the photoemission experiments gave that the effective Hubbard inter- 
action {U) for the Cu ions is of the order U 6 — 10 eV [25], which is much larger 
than the observed band width W (~ 2 eV) [26]. (iii) A rather direct evidence for strong 
correlations comes from the doping dependence of the dynamic conductivity a (ft)) in 
La2-xSrxCu04 and Ndj-xCcxCuO^-y, particularly from the observed shift of the spec- 
tral weight from high to low energies with doping [27]. Besides the development of the 
Drude peak around ft) = in the underdoped systems the so called mid-infrared (MIR) 
peak is also developed around 0.4 eV. 

Regarding the EPI one can put an "old fashioned" question: Does the EPI makes 
(contributes to) the superconducting pairing in HTSC oxides? Surprisingly, most of 
researchers in the field believe that the EPI is irrelevant and that the pairing mecha- 
nism is due to spin fluctuations and strong correlations alone- see [29]. This belief is 
mainly based on an incorrect stability criterion (which, if true, would strongly limit 
Tc in the EPI mechanism), and also on a number of experimental results which give 
evidence for strong anisotropic (d — wave like) pairing with gapless regions on the 
Fermi surface [6], etc. Moreover, the phase sensitive SQUID measurements of the 
Josephson effect [30], [31] in the orthorhombic material YBa2Cu30(,+x are strongly 
in favor of an "orthorhombic" d — wave superconducting order parameter, for instance 
A(k) = As + A(i{coskx — cosky). As experiments of Tsuei et al. [30], [31] show one has 
Ay < O.IA^/ in optimally doped YBa2Cu30(i^x^ which means that zeros of A(k) are near 
intersections of the Fermi surface and the lines kx ~ ±ky. Recent experiments on the 



single-layer crystals Tl2Ba2CuO^+x and on Bi2Sr2CaCu20g^x (Bi22l2) done by Tsuei 
group [32], [33], [34], prove the existence of pure d — wave pairing in underdoped, opti- 
mally and overdoped systems. The recent interference experiments on Nd2-xCexCuO^-y 
point also to d-wave pairing in this compound [35]. In that respect, we point out that 
there is also an widespread (and unfounded) belief that d — wave is incompatible with 
the EPI pairing mechanism. 

Another argument used against the EPI as an origin of superconductivity in HTSC 
oxides is based on the small value of the oxygen isotope effect ao (a = oco + occw + 
OCy -|- asa) in optimally doped materials, such as YBCO with highest critical temperature 
Tc fti 91K where ao ~ 0.05 [105], instead of the canonical value a — 1/2 which would 
be in the case of the EPI pairing mechanism alone and in the presence of O-vibrations 
only. 

On the other hand, there is good experimental evidence that the EPI is sufficiently 
large in order to produce superconductivity in HTSC oxides, i.e. X > I. Let us quote 
some of them: (1) The superconductivity induced phonon renormalization [3], [36], [4], 
[5] is much larger in HTSC oxides than in LTSC superconductors. This is partially due to 
the larger value oi /S./Ef in HTSC than in LTSC; (2) the line-shape in the phonon Raman 
scattering is very asymmetric (Fano line), which points to a substantial interaction of the 
lattice with some quasiparticle (electronic liquid) continuum. For instance, the recent 
phonon Raman measurements [4] on H gBa2Ca2Cu40 lo+x at T < Tc give very large 
softening (self-energy effects) of the Aig phonons with frequencies 240 and 390 cm~^ 
by 6 % and 18 %, respectively. At the same time there is a dramatic increase of the 
line-width immediately below T^, while above T^. the line-shape is strongly asymmetric. 
A substantial phonon renormalization was obtained in (Cu,C)Ba2Ca2Cu40io+x [5]; (3) 
the large isotope coefficients (ao > 0.4) in YBCO away from the optimal doping [105] 
and ao ~ 0.15 — 0.2 in the optimally doped Lai s55'ro.i5CM04. At the same time one 
has ao ~ Ocm making a 0.25 — 0.3. This result tell us that other, besides O, ions 
participate in pairing; (4) the most important evidence that the EPI plays an important 
role in pairing comes from tunnelling spectra in HTSC oxides, where the phonon-related 
features have been clearly seen in the / — V characteristics [37], [38], [39], [40], [41]; 
(4) the penetration depth in the a-b plane of YBCO is increased significantly after the 
substitution ^ o^^, i.e. {Mab/ Xab) = {^^Kb -^^ Kb)I^^Kb = 2.8 %at4K [42]. 
Since X^h ~ fn* the latter result, if confirmed, could be due to the nonadiabatic increase 
of the effective mass m*. 

Recent ARPES measurements on HTSC oxides [43], [109] show a kink in the quasi- 
particle spectrum at characteristic (oxygen) phonon frequencies in the normal and su- 
perconducting state. This is clear evidence that the EPI is strong and involved in pairing. 

On the theoretical side there are self-consistent LDA band-structure calculations 
which (in spite of their shortcomings) give a rather large bare EPI coupling constant 
A ~ 1.5 in Lai 855'ro.i5CM(94 [51], [53]. The nonadiabatic effects due to poor metallic 
screening along the c-axis may increase A additionally [52], [53]. All these facts are 
in favor of the substantial EPI in HTSC oxides. However, if the properties of the nor- 
mal and superconducting state in HTSC oxides are interpreted in terms of the standard 
EPI theory, which holds in LTSC systems, some puzzles arise. One of them is related 
to the normal-state conductivity (resistivity) - in optimally doped systems the width of 
the Drude peak in C7((0) and the temperature dependence of the resistivity p{T) are not 



incompatible with the strong-coupling theory with A ~ 3 and Xtr ~ 1 (if (Opi ~ 3 eV), 
where Xtr is the transport EPI coupling constant [50]. On the other side the combined 
resistivity and low frequency conductivity (Drude part) measurements give ~ 0.3 if 
the plasma frequency takes the value (Qpi ~ 1 - see more below. If one assumes that 
Xtr ^ X, which is the case in most low temperature superconductors (LTSC), such a 
small X can not give large Tc{'^ \QQ K). 

In the past there were doubts on the ability of the EPI to explain the linear temperature 
dependence of the resistivity in the underdoped system [58] Bi2+xSr2-yCuO(,^§, which 
starts at low T > 10 — 20 ^. Because the asymptotic behavior of p (T) (for T <^ ®d) 
is absent in this sample, then it seems that this experiment is questioning seriously the 
contribution of the EPI to the resistivity. However, there are other measurements [59] 
on Bi2+xSr2-yCuOQj.§ where the linear behavior starts at higher temperature, i.e. at 
r > 50 ^. Additionally, the resistivity measurements [14] on Bi2SrCuOx samples with 
low Tc^3 K show saturation to finite value at T = ^. After subtraction of this constant 
part one obtains the Bloch-Griineisen behavior between Tc^3 K and 300 K, which is 
due to the EPI. 

Concerning the EPI, the above results imply the following possibilities: (a) A^,. ^ 1 < 
X and the pairing is due to the EPI, or (b) Xtr — X^ 0.4 — 0.6 and the EPI is ineffective 
(although present) in pairing; (c) Xtr — X but the EPI is responsible for pairing on the 
expense of some peculiarities of equations describing superconductivity. In Section 5. 
we present a theory of the EPI renormalized by strong electronic correlations, which is 
in favor of the case (a). It is interesting that the similar puzzling situation (Xtr ^ X) 
is realized in BaxKi^xBiO^, compound (with Tc ~ 30 K), where optical measurements 
give Xtr ~ 0.1 — 0.3 [15], while tunnelling measurements [16] give A ~ 1. Note, in 
BuxKi^xBiO^ there are no magnetic fluctuations (or magnetic order) and no signs of 
strong electronic correlations. Therefore, the EPI is favored as the pairing mechanism 
in BaxKi-xBiO^. It seems that in this compound long-range forces, in conjunction with 
some nesting effects, may be responsible for this discrepancy? 

One can summarize, that the EPI theory, which pretends to explain the normal metallic 
state and superconductivity in HTSC oxides, is confronted with the problem of explain- 
ing why the EPI coupling is present in self-energy effects (governed by the coupling 
constant A > 1) but it is suppressed in transport properties (which depend on Xtr < 1), 
i.e. why A^r is (much) smaller than A. One of the possibilities is that strong electronic 
correlations, as well as the long-range Madelung forces, affect the EPI significantly. This 
will be discussed in forthcoming sections. In light of the above discussion it is also im- 
portant to know the role of the EPI in the formation of d — wave superconducting state 
in HTSC oxides, i.e. why it is compatible with d-wave pairing? 

In this review we discuss theoretical and experimental results in HTSC oxides and 
mostly those which are related to: (i) strong quasiparticle scattering in the normal state, 
(ii) the pairing mechanism [17], [18], [19], [20]. 

The paper is organized as follows. In Section 2. we review important physical proper- 
ties of HTSC oxides in the normal and superconducting state, whose understanding is a 
basis for the microscopic theory of superconductivity. Only those experiments (and the- 
oretical interpretations) are discussed here which are in our opinion most important in 
getting information on the pairing mechanism in HTSC oxides. In Section 3. we discuss 



the general theory of the EPI and its low-energy version. The theory of strong electronic 
correlations is studied in Section 4., where much space is devoted to a systematic, re- 
cently elaborated, method for strongly correlated electrons [17], [18], [19], [20] - the 
X-method. The latter considers strongly interacting quasiparticles as composite objects, 
contrary to the slave-boson method which at some stage assumes spin and charge sep- 
aration [28]. A systematic theory of the renormalization of the EPI coupling by strong 
electronic correlations [17], [18] is exposed in Section 5. It is shown there, that tha for- 
ward scattering peak develops in the EPI by lowering doping, while the coupling at 
large transfer momenta (the backward scattering) is suppressed . 

In Section 6. we summarize the basic predictions of the theory based on the existence 
of the forward scattering peak in the EPI, impurity and Coulomb scattering, and possible 
relation between the forward scattering peak in the EPI and pseudogap. The comparison 
between the EPI and SFI prediction is given in Section 7. The (im)possibility of super- 
conductivity in the Hubbard and t-J mode is studied in Section 8., while the obtained 
results are summarized in Section 9. 



1.2. Prejudices on the EPI 

In spite of the reach experimental evidence in favor of the strong EPI in HTSC 
oxides there was a disproportion in the research (especially theoretical) activity, since 
the investigation of the spin fluctuations mechanism of pairing prevailed in the literature. 
This was partly due to a theoretically unfounded statement - given in [10], on the upper 
limit of Tc in the phonon mechanism of pairing. It is well known that in an electron- 
ion system besides the EPI there is also the repulsive Coulomb interaction and these 
are not independent. In the case of an isotropic and homogeneous system with a weak 
(quasi)particle interaction the effective potential Ve//(k, (o) in the leading approximation 
looks like as for two external charges (e) embedded in the medium with the total 
macroscopic longitudinal dielectric function e;o;(k, Ct)) (k is the momentum and (O is 
the frequency) [12], i.e. 

In the case when the interaction between quasiparticles is strong, the state of embed- 
ded quasiparticles changes substantially due to the interaction with other quasiparti- 
cles, giving rise to Ve//(k, (O) ^ Ane^ /k^Etoti}^, (o). In that case Veff depends on other 
(than er„f(k, O))) response functions. However, in the case when Eq.(l) holds the weak- 
coupling limit is realized where T^ is given by Tc — ft)exp(— 1/(A — jU*) [9], [12]). 
Here, A is the EPI coupling constant, (b is the average phonon frequency and ju* is the 
Coulomb pseudo-potential, n* = / {I + n InEp/a) (Ep is the Fermi energy). A and jU 
are expressed by £tot{K co = 0) 

rlkp hf]h Attp^ 



where A'(O) is the density of states at the Fermi surface and kp is the Fermi momentum - 
see more in [1 1]. In [10] it was claimed that the lattice stability of the system with respect 
to the charge density wave formation implies that the condition etot (k, ft) = 0) > 1 must 
be fulfilled for all k. If this were correct then from Eq.(2) follows that }l > X, which 
limits the maximal value of Tc to the value T^^^ ~ Ep exp(— 4 — 3/ A) . In typical metals 
£'f < (1 — 10) eV and if one accepts this (unfounded) statement that A < jU < 0.5 one 
obtains ~ ( 1 — 10) ^. The latter result, of course if it would be true, means mean that 
the EPI is ineffective in producing high-T^ superconductivity, let say not higher than 20 
Kl However, this result is apparently in conflict with a number of experimental results 
in low-Tc superconductors (LTS), where jU < A and X > I. For instance, X a; 2.5 is 
realized in PbBi alloy, which is definitely much higher than /i(< 1) thus contradicting 
the statement made in Ref.[10]. 

The statement in [10] that eun (k, £0 = 0) > 1 must be fulfilled for all k is in an apparent 
conflict with the basic theory [12], which tells us that £fof(k, (o) is not the response 
function. If a small external potential 5Vext{^-,0i) is applied to the system it induces 
screening by charges of the medium and the total potential is given by SV^of (k, to) = 
5Vex?(k, (o)/e/„f(k, (o) which means that l/gfo? (k, (o) is the response function. The latter 
obeys the Kramers-Kronig dispersion relation which implies the following stability 
condition: l/ejof (k, O) = 0) < 1 for k 7^ 0, i.e. either efor(k 7^ 0, w = 0) > 1 or Cfof (k ^ 
0, <o = 0) < 0. This important theorem has been first proved in the seminal article by 
David Abramovich Kirzhnits [12] and it invalidates the formula for by setting 
aside the above restriction on the maximal value of T^. 

Is efo/(k 7^ 0,(0 = 0) < realized in real systems? This question was thoroughly 
studied in Ref. [13] and in the context of HTSC in [11], while here we enumerate 
the main results. In the inhomogeneous system, such as a crystal, the total longitu- 
dinal dielectric function is matrix in the space of reciprocal lattice vectors (Q), i.e. 
ef(,f(k + Q,k-|-Q',ft)), and ero?(k, w) is defined by £jo/(k, w) = £fo/(k + 0,k + 0, ft)). 
For instance in the dense metallic systems with one ion per cell (such as the metallic 
hydrogen ) and with the electronic dielectric function eg/(k, 0) one has [13] 

e,,,(k,0) = — ^^^^1^. (3) 

e,,(k,0)G£p(k) 

At the same time the frequency of the longitudinal phonon ft)/ (k) is given by 

""'^^^ = ^^^^ ^""^ 

where Gep is the local field correction Gep - see Ref. [13]. The right condition for the 
lattice stability requires that the phonon frequency must be positive, ft)f (k) > 0, which 
implies that for ee/(k,0) > onehas £e/(k,0)G£p(k) < l.The latter gives £fof(k,0) < 0. 
The calculations [13] show that in the metallic hydrogen crystal £fof(k, 0) < for all 
k 7^ 0. The sign of £fo;(k, 0) for a number of crystals with more ions per unit cell is 
thoroughly analyzed in [13], where it is shown that e?o<(k 7^ 0,0) < is more a rule 
than an exception. The physical reason for £;«? (k 7^ 0,0) < is due to local field effects 
described by GEpiS). Whenever the local electric field E/^c acting on electrons (and 



ions) is different from the average electric field E, i.e. E/^c 7^ E there are corrections 
to ero,(k,0) (or in the case of the electronic subsystem to ee(k,0)) which may lead to 
etotiKO)<0. 

The above analysis tells us that in real crystals e/o;(k,0) can be negative in the large 
portion of the Brillouin zone giving rise to A — /i > 0, due to local field effects. This 
means that the dielectric function Etot does not limit Tc in the phonon mechanism of 
pairing. The latter does not mean that there is no limit on Tc at all. We mention in 
advance that the local field effects play important role in HTSC oxides, due to their 
layered structure with ionic-metallic binding, thus giving rise to large EPI - see more 
subsequent sections. 

In concluding we point out, that there are no theoretical and experimental arguments 
for ignoring the EPI in HTSC oxides. However, it is necessary to answer several impor- 
tant questions which are also related to experimental findings in HTSC oxides: (1) If the 
EPI is responsible for pairing in HTSC oxides and if superconductivity is of J — wave 
type, how these two facts are compatible? (2) Why is the transport EPI coupling constant 
Xtr (entering the resistivity formula) much smaller than the pairing EPI coupling con- 
stant A (> 1 ) (entering the formula for Tc), i.e. why one has 0.4 - 0.9) < A (~ 2)? 

(3) Is the high Tc value possible for a moderate EPI coupling constant, let say for A < 1? 

(4) Finally, if the EPI interaction is ineffective for pairing in HTSC oxides why it is so? 



2. EXPERIMENTS RELATED TO PAIRING MECHANISM 

A much more extensive discussion (than here) of the experimental situation in HTSC 
oxides is given in a number of papers - see reviews [2], [1 1] . In the following we discuss 
briefly experimental results, by including the most recent ones, which can give us a clue 
for the pairing mechanism in the HTSC oxides. 



2.1. Magnetic neutron scattering 

2.7.7. Normal state 

The cross-section for the inelastic neutron magnetic scattering is expressed via the 
Fourier transform of the spin-correlation function (the spin structure factor) 5""(k, ca) 
which is proportional to the imaginary part of the susceptibility Imx{Kk^,0)). In the 
(normal) metallic state of doped HTSC oxides without magnetic order the inelastic 
scattering (in absence of the AF magnetic order) is of interest and in most systems 
/m;^(k, is peaked around the AF wave- vector Q = {n, k). The pronounced mag- 
netic fluctuations in the underdoped metallic state is contrary to usual metals (described 
by the Landau-Fermi liquid) where the magnetic fluctuations are much weaker. In HTSC 
oxides /m;t;(k, (o) depends on hole doping, and for instance, in the bilayer (two layers 
per the unit cell) compound YBa2Cu-iO(,^^ the low energy spectra is peaked at Q, whose 
width 5m broadens by increasing doping concentration - see review [66]. Around the 
optimal doping the magnetic correlation length = (2/ 5m) ~ (1 — 2)a is almost tem- 



perature independent. This fact contradicts the assumption of the theory of spin fluctu- 
ation mechanism by the Pines group [29], where is strongly T-dependent. We stress 
that in the SFI theory Imx (k, ft)) is important quantity since the effective pairing po- 
tential Veff(li., (o) and the self-energy Lj/(k, a) are approximately given by (on the real 
frequency axis - see [2]) 

—G{q,a)Veff{k + q,(0 + Q), 
2 daimxiq,^ + iO^) 



VsF{q,(0 + iO+) = 8lF — ^ (5) 

J oo n 12 — ft) 

where G(q,^2) is the electron Green's function. This approach can be theoretically 
justified in the weak coupling limit (U <^ W) only . Although the HTSC oxides are 
far from this limit this expression is frequently used in the SFI theories of pairing, where 
larger /m;^(k, ft)) should give larger T^. 

What is the experimental situation? The antibonding (odd) spectral function 
Im.x^"'^^\k,(o) of YBa2Cu^0(,+x is strongly doping dependent as it is seen in Fig. 1. 
By comparing the magnetic neutron scattering (normal state) spectra in F5a2CM3 06.92 
and YBaiCu^O in Fig. la the difference is reflected in their spectral functions 
Imx^"^'^'^ (k, ft)) . Namely, in the frequency interval which is important for superconduct- 
ing pairing Imx^°''^^ (Q, co) of 

YBa2Cu20e.92 is much larger than that in YBa2CuzO(,_g'] although the differences 
in their critical temperatures T^- is very small, i.e. Tc = 9\ K for YBa2Cu2,0^ i)2 and 
Tc = 92.5 K for YBa2Cu^0(yS)i- This result, in conjunction with the anti-correlation 
between the NMR spectral function /q = \mi(i)^Qlmx{Q,(o) / (O and T^ - shown in 
Fig. la, is apparently against the SFI theoretical models for pairing mechanism [29], 
[65]. 



2.7.2. Superconducting state 

In the superconducting state the magnetic fluctuations are drastically changed, what 
is in fact expected for the singlet pairing state which induces spin gap in the magnetic 
excitation spectrum of s-wave superconductors. However, the spectrum in the supercon- 
ducting state of HTSC oxides is more complex due to d-wave pairing and specificity 
of the band structure. For instance, at T < Tc the sharp peak in Imx^"'^'^HK ft)) is seen 
at COreson = 41 meV and at k2D = (n/a^n/a) of the fully oxygenated (optimally doped) 
YBa2Cu20^+x (-^ ~ 1; pa 92 K) [71], [72]. The doping dependence of the peak posi- 
tion and its width [66] is shown in Fig. Ic, where it is seen that by increasing doping 
the peak in the superconducting state becomes sharper and moves to higher frequencies 
(scaling with T^), while its height is decreasing. This can be qualitatively explained by 
using the RPA susceptibility 




t>/> Emrgy (meV) 

FIGURE 1. Magnetic spectral function Imx'<~'i (k, ©): (a) Iq{Tc) values at T = 200 K for various HTSC 

oxides: LSCO - La2xSr^CuO^\ TBCO - ThBaiCuOe+x and Tl2Ba2CaCu20g; YBCO - and YBa2Cu40g - 
from [67]; (b) for YBa2Cuj,0^+x in the normal state at T = 100 K and at g = 100 counts in the 

vertical scale corresponds to Xmwc ~ SSOjuJ/eV - from [66]; (c) for YBa2CuT,0(,+x in the superconducting 
state at r = 5 A" and at 2 = (;r, ;r) - from [66]. 



where the bare susceptibiUty ;t;o(k, co) contains the coherence factor [1 — ((^k+q<^q + 
Ak+qAq)/£'k_|_q£'q] - see [2]. This (type II) coherence factor reflects the (well known) fact 
that the magnetic scattering is not the time reversal symmetry. In the case when k and k + 
q are near the Fermi surface and when A^+q pa — Aq at k = Q = (;r/(5[, :;r/a) the coherence 
factor is of the order of one at or near the Fermi surface (note (^k+q<^q < 0) and therefore 
contributes significantly to ;to(k = Q, (o). The case AqAk+q < is realized when the 
d — wave order parameter, for instance A^ = ( Ao/2) [cos — cos . So, the mechanism 
of the peak formation (below T^) is the consequence of the electron-pair creation with 
an electron in the (+ ) lobe and a hole in the (— ) lobe of the superconducting order 
parameter. Note, that the (±) lobes of A^ are separated approximately by the wave- 
vector Q = (^r/a, k/o). Due to the large density of states near the lobes a large peak in 
Imxi^ — Q) Ki is expected to be realized, i.e. COreson > 2Ao. Of course the better (than 
RPA) calculations of ;^(k = Q,^^, w) is needed for a full quantitative analysis, where 
a possible resonance in x{K(o) with COreson < 2Ao can also contribute. It is important 
to stress, that the magnetic resonance in the superconducting state is consequence of 
superconductivity but not its cause as it was stated in some papers. It can not be the 
cause for superconductivity simply because its intensity at T around T^. is vanishing 
small and not affecting Tc at all. If the magnetic resonance would be the origin for 



superconductivity (and high T^) the phase transition at must be first order, contrary 
to experiments where it is second order. 

The next very serious argument against the SFI pairing mechanism is the smallness of 
the coupling constant ggf. Namely, the real spin-fluctuation coupling constant is rather 

small gsf < 0.2 eV, what is in contrast to the large value (g^J^^^^ ~ 0.6 eV) assumed 
in the SFI theory by the Pines group - see the MMP model in Section 7.1. The upper 
limit of gsf{< 0.2 eV) is extracted from : (i) the width of the resonance peak [68], and 
(ii) the small magnetic moment (jU < 0. 1 jUg) in the antiferromagnetic state of LASCO 
and YBCO [69]. Note, that the pairing coupling in the SFI theory is A,/ ~ g^^, and for 
the realistic value of gsf < 0.2 eV it would produce A^y ~ 0.2 and very small Tc ^ I K. 
The SFI model roots on its basic t — J Hamiltonian. However, recently it was shown in 
[70] that there is no superconductivity in the t-J model at temperatures characteristic for 
HTSC oxides - see Fig. 27 below. If it exists T^. must be very low. 

In conclusion, the inelastic magnetic neutron scattering give evidence that the spin 
fluctuations interaction (SFI), although pronounced in underdoped systems, is inejfective 
in the pairing mechanism of HTSC oxide. However, the SFI in conjunction with the 
residual Coulomb repulsion triggers superconductivity from s-wave to d-wave, whose 
strength is predominantly due to the EPI - see discussion in Sections 5.-7.. 



2.2. Dynamical conductivity and resistivity p{T) 

Since a{G)) and p{T) give important information on the dominant scattering mecha- 
nism, in the following we analyze their properties in more details. 



2.2.1. Dynamical conductivity <y{(o) 

o{co) is in fact derived quantity since it is extracted from the measured optic reflectiv- 
ity R{(0) and absorption A (co). By measuring the normal-incident (of light) reflectivity 
R{(o) in the whole frequency region (0 < (0< <») one can determine the phase {(o) of 
the complex reflectivity 

r(a)) = = (7) 

Ve(<o) + 1 

by the Kramers-Kronig relation, and accordingly to determine in principle the complex 
dielectric function 

/ N 4;ri(7((o) 

= eoo + e/a«M + — (8) 

where Ccx, and (7((u) are electronic contributions and eian is the lattice contribution. 
However, R{(o) is usually measured in a finite (O region and extrapolations is needed, 
especially at very low frequencies. This extrapolation of R{(o) also contains some model 
assumptions on the scattering processes in the system (on a{G))), i.e. 1 —R{(o) ~ 
- the Hagen-Rubens relation for the standard (with elastic scattering only) Drude metal. 



or 1 — ~ ft) for strong EPI (or for marginal Fermi liquid). So, one should be always 
cautious not to overinterpret the meaning of a (ft)) obtained in such a way. 

In HTSC oxides R{(o), A{(o) are usually measured in a broad frequency region - up 
to several eV. At such high frequencies the interband transitions take place and in order 
to calculate o{0)) the knowledge of the band structure is needed. This problem was 
analyzed in the framework of the LDA band structure calculations [73] by taking into 
account the interband transitions, where a rather good agreement with experiments for 
0) > I eV was found. This is surprising since the LDA-method does not contain the 
Hubbard bands, and according to [73] there is no sign of transitions between Hubbard 
sub-bands in the high energy region of a (ft)). This very interesting result deserves to be 
further analyzed since it contradicts the physics of the Hubbard models. 

Here we discuss briefly the normal state a (ft)) in the low frequency region (O < I eV 
where the intraband effects dominate the quasiparticle scattering. In the low ft) regime 
the processing of the data in the metallic state of HTSC oxides is usually done by using 
the generalized Drude formula for the inplane conductivity a (ft)) = <7i + i02 [86], [87], 
[88] 

(^l a 1 
~ An Ttr{(0j)-i(0m,r{(0)/mj 

i = a,b enumerates the plane axis, Ttrico.T) and m„-{G)) are the transport scattering 
rate and optic mass, respectively. Sometimes in the analysis of experimental data the 
effective transport scattering rate T*f^{(0,T) and the effective plasma frequency ft)p(ft)) 
are used, which are defined by 

mtr{CO) O2{C0) 



and 

<(ft)) = ^— ft)2. (11) 

mtr{(0) ' 



For best optimally doped HTSC systems the best fit for rf^(G), T) is given by Tf^{0), T) 
max{ar,j8ft)} in the temperature and frequency range from very low (~ 100 K) up to 
2000 K, where a,j8 are of the order one - see Fig. 2. These results tell us that the 
quasiparticle liquid, which is responsible for transport properties in HTSC, is not a sim- 
ple (weakly interacting) Fermi liquid. We remind the reader that in the usual (canoni- 
cal) normal Fermi liquid with the Coulomb interaction on has rf^(ft), T) ~ ^^(ft), T) ~ 
r(ft),r) ~ max{r^, O)^} at low T and ft), which means that quasiparticles are well de- 
fined objects near (and at) the Fermi surface since ft) ^ r(ft), T). In case of HTSC oxides 
with r(ft), T) ~ max(ft), T) in the broad regions and the quasiparticles decay rapidly and 
therefore are not well defined objects. At these temperatures and frequencies the simple 
canonical Landau quasiparticle concept fails. The latter behavior can be due to the strong 
electron-electron inelastic scattering, or due to the quasiparticle scattering on phonons 
(or on other bosonic excitations). It is important to stress, that quasiparticles 

interacting with phonons at finite T are not described with the standard Fermi liquid, 
in particular at T > 0z)/5, since the scattering rate is larger than the quasiparticle 
energy, i.e. one has F ~ max(ft),r). Such a system is well described by the Migdal- 
Eliashberg theory whenever (Od <C Ep is fulfilled, which in fact treats quasiparticles 
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FIGURE 2. The transport scattering rate \/t{(0) (in the text rtr(ft))) and the transport effective mass 
m*{(o)/me (in the text mtr{(£>) / moo) for series of underdoped HTSC oxides. Ttr{(i>) is temperature inde- 
pendent above 1000 cm^^ but it is depressed at low T and low (o- from [62]. 



beyond the original Landau quasiparticle concept. Note, that even when the original 
Landau quasiparticle concept fails the transport properties may be described by the 
Boltzmann equation, which is a wider definition of the Landau-Fermi liquid. 

We point out, that in a number of articles it was incorrectly assumed that r((0, T) 
Ttr{(O^T) ^ Y^^{(o,T) holds in HTSC oxides . The above discussed experiments (see 
Fig. 2) give that Tj^{(0,T) is linear in the broad region of CO and T up to 2500 K - see 
in Fig. 2. However, if Y{(O^T) is due to the EPI it saturates at the maximum phonon 

frequencies (BSax(< 1000 K). By assuming also that T^"{(0, T) ^ r^^^((0, T) holds for 
all CO, in a number of papers it was concluded, that the EPI does not contribute to the 
inelastic scattering of quasiparticles and to the Cooper pairing in HTSC oxides. Does it 
hold r^^'{(0, T) fti rf/^(a), T) in HTSC oxides? The answer is NO. 

o{(o) of HTSC oxides was theoretically analyzed [87], [88] in terms of the EPI, 
where it was found that Ttr{(0,T) and mtr{0),T) depend on the transport spectral 
function af^^pF{(o) - see more in [2]. Their analysis is based on: (i) the assump- 
tion that af^^pF{(o) ~ a^pF{co) - the Eliashberg spectral function: (ii) the shape of 

a^pF{(o) is extracted from various tunnelling conductivity measurements, [38], [39], 
[40], [41], which makes a rather large EPI coupling constant and the critical temperature 
X — 2 d(oa?'{(o)F{(i)) /co ^ 2 and Tc ^ 90 K, respectively; (iii) the plasma frequency 
is taken to be COpi = 3 eV. It was obtained that Tf/'^O), T) ~ Ct) in a very broad ft)-interval 

(up to 250 meV), which is much larger than the maximum phonon frequency £o£ax ~ 80 
mev. This is illustrated in Fig. 3. Moreover, rf/'{o),T) differs significantly from the 
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FIGURE 3. The theoretical predictions for the frequency dependence of the various relaxation rates 
y(= r) with a^F{(o) - right: the generalized Drude fit for rtr{(o) - solid line; rj'r(co) - short-dashed fine; 
r{(o) - long dashed line; r,r(aj) calculated - dotted line. In the inset the calculated r;*^(aj)(= 1/t*{(o)) 
with (sohd line) and without (dashed line) the interband contributions with a^F{(0) from right and at 
r= 100/: -from [88]. 

quasiparticle scattering rate r^^((0, T) = — 2/mE(ft)) [87], [88]. We see from Fig. 3 that 
r^^((U,r) is much steeper function than rf/'{o),T) and the former saturates at much 
lower frequency - of the order of the maximum phonon frequency C!)max- 

Note, that r*f {(O, T) = {m/mtr{(o))Tff{(o, T) is also quasi-linear function in a very 
broad region \5QK<(0< 3000 K - see Fig. 3 . The slope of T*^^^ ( ct) , T ) is of the order of 
one, in accordance with experiments results [87], [88], and it (and r^^((0, T)) saturates 

at (Osat — —ImLtr{(Osat) ^ tt^Sax Only. The transport spectral function a^^{(o)F{(o) can 
be also extracted from the transport scattering rate rtr{o},T = 0) - see [2], [11], since 
the theory gives that 



However, real measurements are performed at finiter(> Tc) where a},,{G))F {o)) is the 
solution of the Fredholm integral equation (of the first kind) . Such an inverse problem 
at finite temperatures in HTSC oxides is studied first in [87] (see also [88]), where the 
smeared structure of af^{(o)F{(o) in YBa2Cu20j-x was obtained, which is in qualitative 
agreement with the shape of the phonon density of states F(ft)) . At finite T the problem is 
more complex because the fine structure of a^j.{(o)F{(D) gets blurred as the calculations 
in [74] show. The latter gave that a^y{co)F{co) ends up at cOmax ~ 70 — 80 meV, which is 
the maximal phonon frequency in HTSC oxides. This result indicates strongly that the 
EPl in HTSC oxides is dominant in the IR optics. We point out, that if R{(o) (and o{(0)) 
are due to some bosonic process with large frequency cutoff (Oc in the spectrum, as it is 




(12) 



the case with the spin-fluctuation (SFI) scattering where COc ~ 400 meV, the extracted 
OCf^(tt))F((o) should end up at this high ©c- The latter is not seen in optic measurements 
at r > Tc, which tells us that the SFI scattering, with a1^{co)F{co) ~ g]j-ImXs{o}) and 
with the cutoff (Oc > 400 meV, is rather weak and ineffective in optics of HTSC oxides.. 

We stress that the extraction of F^^ from R{(0) is subtle procedure, because it depends 
also on the assumed value of e^o. For instance, if one takes Coo = 1 then Tf/ is linear 
up to very high CO, while for Coo > 1 the linearity of Tf/' saturates at lower a. Since 
rf/'{(D, T), extracted in [62], and recently also in [63], is linear up to very high (O it may 
be that the ion background and interband transitions (contained in e,^) are not properly 
taken into account in these papers. As a curiosity in a number of papers, even in the 
very cited ones such as [62], [63], there is no information which value for Coo they take. 
We stress again, that the behavior of Ttr{(o) is linear up to much higher frequencies for 
Coo = 1 than for Ecx, ~ 4 — 5 - the characteristic value for HTSC, giving a lot of room 
for inadequate interpretations of results. In that respect, the recent elipsometric optic 
measurements on YBCO [75] confirm the results of the previous ones [76] that £00 > 4 
and that TfJ^ saturates at lower frequency than it was the case in Ref. [62]. We stress 
again that the reliable estimation of the value and (O, T dependence ofYtr{o)) and m{o)) 
can be done, not from the reflectivity measurements [62], [63], but from elipsometric 
ones only [76], [75]. 

In concluding this part we stress two facts: (1) The large difference in the (O^T 
behavior of Tff{o},T) and T^^{(o,T) is not a specificity of HTSC oxides but it is 
realized also in a number of LTSC materials. In fact this is a common behavior even 
in simple metals, such as Al, Pb, as shown in [118], where T^^{(0,T) saturates at 
much lower (Debay) frequency than Tf^{(O^T) and T*f'^{(o,T) do. In that respect 
the difference between simple metals and HTSC oxides is in the scale of phonon 

frequencies, i.e. (O^^x ~ 100 K in simple metals, while £0(^ax ~ 1000 K in HTSC oxides. 
Having in mind these well established and well understood facts, it is very surprising 
that even nowadays, 18 years after the discovery of HTSC oxides, the principal and 
quantitative difference between F and Ttr is neglected in the analysis of experimental 
data. For instance, by neglecting the pronounced (qualitative and quantitative) difference 
between Ttr{0), T) and r((0, T), in the recent papers [63], [64] were made far reaching, 
but unjustified, conclusions that the magnetic pairing mechanism prevails; (2) It is worth 
of mentioning, that quite similar (to HTSC oxides) properties, of o{(o), R{(o) and 
p(r) were observed in experiments [92] on isotropic metallic oxides Lao ^Sro sCoOs 
and Cao ^SrQ sRuOs - see Fig. 4. We stress that in these compounds there are no signs 
of antiferromagnetic fluctuations (which are present in HTSC oxides) and the peculiar 
behavior is probably due to the EPI. 

It is worth of mentioning that after the discovery of HTSC in 1986 a number of contro- 
versial results related to o{(0) were published, followed by a broad spectrum of results 
and interpretations, from standard approaches up to highly exotic ones. For example, the 
reported experimental values for (Opi were in the surprisingly large range (0.06 — 25) eV, 
causing a number of exotic (and confusing) theoretical models for electronic dynamics 
- see more in [76]. (The similar situation was with ARPES measurements - see below.) 
So, one should be very cautious in interpreting experimental and theoretical results. In 
that respect, recent experiments related to the optical sum-rule is an additional example 
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FIGURE 4. Broad range specular reflectance spectra of CaosSrosRuO^ (broken line) and 
LaosSro^^CoO^ (solid line). Inset spectra of Tl2Ba2Ca2Cuj,0\Q, Bi2Sr2CaCu20s, YBa2Cu^0-] and 
La\,^sSro.\5CuOA. From [92]. 



for controversies in this field coming from inadequate interpretations of results. This is 
the reason why we devote more space to the problem of "violation" of partial sum-rule. 

There are two kinds of sum rules which are used in interpreting results on a (ft)). The 
first one is the total sum rule and in the normal state it reads 

/ af G) Jo) = ^ = -— , (13) 

,/o 8 2m 

while in the superconducting state [77] it is given by the Tinkham-Ferrell-Glover (TFG) 
sum-rule 

-^^j^aU<o),. = -l-. (14) 

Here, n - is the total electron density, e - the electron charge, m - the bare electron 
mass, Al - the London penetration depth. The first term c^/8A^ is due to the appearance 
of the superconducting condensate (ideal conductivity) which contributes <7j {co) = 



(c^/4A^)5(tt)). The total sum rule represents the fundamental property of matter - the 
conservation of the electron number. To calculate it one should use the total Hamiltonian 
Htot = Te + Hint by taking into account all electrons, bands and their interactions Hint 
(Coulomb, EPI, with impurities,etc.). Here Te is the kinetic energy of bare electrons 

fe = Y.j d\\ifl{x)^\if^{x) = £ ^4o^po. (15) 

The partial sum rule is related to the energetics in the conduction (valence) band. 
Usually it is derived by using the Hamiltonian of the valence electrons 

Hv = Tv + Vy^coul = 2^ ^p(^v,pa^v,pa + ^vfiouU (16) 

p,C7 

which contains the band-energy (with dispersion gp) and the Coulomb interaction of 
valence electrons Yv^Coul- In the normal state the partial sum-rule reads [78] (for general 
form of Cp) 

■<«>,M» = |^£<!^.™^ (17) 



^0 



where nv,p = Cpcr^pa and the reciprocal mass is given by 1 /mp = d^e^j dp\.To simplify 
further discussion we assume for gp = —lt{co%pxa^co^pya) the tight-binding model 
with nearest neighbors (n.n.) where 1/mp = — 2ra^cos Pxa. In practice measurements are 
performed up to finite ft) and the integration over ft) goes up to some cutoff frequency (Oc 
(of the order of band plasma frequency). It is straightforward to show that one has (for 
the n.n. tight-binding model) 



L 



. <v(<oM<» = -^{-Tv) = (18) 

2 o 



where {—Tv)h,. = — Lp£p(«v)/fi, and by ft)^; ^ is defined the band plasma frequency. 
In that case the partial sum-rule in the superconducting state reads 



^ + cTt(a))rfa) = -^(-r,). (19) 



The sum-rule was studied intensively in optimally and underdoped Bi^Sr^CaCu^O^^x 
and in YBa2Cu^0q x for the intraplane conductivity, where the whole frequency region 
is separated into the low ("intraband")- and high ("interband")-frequency parts Al and 
{Ah +Avh), respectively 

Al(0, ft)c) +A//(ft)c, acoc) -\-Avh{oc(Oc,°°) = (Opi/S (20) 

with Al(0, ft)t ) = Ai(0, ft)t ) + ^SNfO^i s where 5sN = 1 in superconducting state and 
A(ft)i,ft)2) = /^ai(ft))Jft). The temperature dependence of Ai and Ah in the above 
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FIGURE 5. Measured T-dependenceof Al(0, (Oc) and Ah {o)c,2o)c) for (Oc « 1.25eV oiBillll {T^ = 
/sT). The data from [79] 



HTSC oxides was studied in [79] and [75], by assuming that "intraband" effects are 
exhausted for ~ 1 - 25 and the main temperature dependence of the high-frequency 
region comes for a = 2 in Eq.(22), i.e. for 1.25 eV < ft) < 2.5 eV , while the temperature 
dependence of the very high energy part Ay// is negligible. It was found that A/^(0, ftJt ) 
grows (quadratically with T) about 6% between 300 and 4 K, while Ah{o3c,occOc) 
decreases with decreasing T^. In the superconducting state there is a small extra increase 
of Al(0, (Dc). The results are shown in Fig. 5. 

In connection with this experiment let us stress that in the BCS superconductor the 



TFG sum-rule is practically satisfied if the integration goes up to (Os ~ (4 — 6)A, where 
A is the superconducting gap. This means that in the BCS superconductors the spectral 
weight appearing in the condensate (at 0) = 0) is transferred from the region 0+ — (B[j. 
However, the experiment in [79] shows a transfer of the spectral-weight from the high 
(£0 > 1 eV) to low energies - see below. This fact was interpreted by some researchers 
[80], [81] as a "violation" of the TFG sum-rule, i.e. that there is more spectral weight 
in the condensate (at ft) = 0) than it is expected from the TFG sum-rule and effectively 
means the decrease of the kinetic energy in the superconducting state. This is in contrast 
to the increase of the kinetic energy in the BCS superconducting state. We are going to 
discuss this problem in details and to demonstrate that the analysis in terms of the kinetic 
energy only, is untenable. 

What is the origin of the spectral-weight transfer, especially in the superconducting 
state of HTSC oxides? Here, we shall study the inplane Oa-b{0)) only, since the origin 
of the quasiparticle dynamics along the c-axis is still unclear. 

The first theoretical interpretation of the spectral-weight transfer was based on the 
partial sum-rule in which the temperature dependence is related to the temperature 
change of the kinetic energy (—71,) (or for more realistic spectrum of co^^ yiT)) - see 

Eq. (17-19). In this framework the extra increase of A/^(0, ft)c) is related to the lower- 
ing of the band kinetic energy in the superconducting state [80], [81]. If this would 
be true, then the lowering of the band kinetic energy (per particle) is approximately 
{{Tv)n,t>Tc -- {Tv)s,T<Tc)/N ~ 1 meV, what is approximately by factor ten larger than 
the superconducting condensation energy. Note, that in the weak coupling BCS theory 
of superconductivity the kinetic energy is increased in the superconducting state. So, the 
alleged large lowering of the kinetic energy in the superconducting state is interpreted 
as a result of some exotic pairing mechanism in which the kinetic energy (or ftj^^ ^{T) 
for more general spectrum) is significantly lowered but the potential energy is increased 
in the superconducting state, contrary to the case of BCS approach. However, this inter- 
pretation misses a very important contribution to the partial sum rule, which is due to 
the large and strongly T, CO dependent transport scattering rate TtriT, (o). 

Before discussing the partial sum rule more adequately, let us mention that the sep- 
aration of the valence-band kinetic energy from the potential one in strongly correlated 
systems is not well defined procedure. For instance, in the Hubbard model with U »t 
and nearest neighbor hoping t one has (see below and also in [2]) the low-energy (va- 
lence) Hamiltonian Hy is given 

H, = -tY, XrX+S^ (21) 
i,5,a 

where the Hubbard operators Z; describe the motion of composite quasiparticles with 
excluded doubly occupancy - see more in Section 4. They have complicated non- 
canonical (anti)commutation rules, which means that Eq.(21) mixes the kinetic energy 
with the (kinematical) potential energy of band (valence) quasiparticles. 

The second theoretical approach is proposed recently [82] , which is in principle exact, 
is based on the fact that in HTSC oxides there is strong electron scattering - direct or 
via phonons, on impurities, etc. So, the presence of an inelastic (and elastic) scattering 
prevents the interpretation of the partial sum rule in terms of the band kinetic energy 



only. As an illustrative example for this assertion may serve the scattering of electrons 
on impurities, where the intraband contribution to Oi^y{(0) is given by 



<yi,v{co) 



47t +r? 



(22) 



Here, Tt^tr/^ — ^hi,tr is the quasiparticle transport relaxation rate due to impurities. In 
this case the partial sum-rule reads 



This result means that the intraband sum-rule can be satisfied in the presence of impuri- 
ties only for ©c oo. The similar conclusion holds in the case of inelastic scattering via 
phonons although in that case Tfr is (O- and T-dependent. The similar reasoning holds for 
interband transitions which in the presence of scattering have also the low-frequency tail. 
Since in HTSC oxides Ytr{0)) is dominantly due to the EPI and reaches values up to 100 
meV , there is no other way to study the partial sum-rule (the value of O\^v{(o)d(o) 
than to calculate O\^v{(o) directly from a microscopic model. The latter must incorporate 
relevant scattering mechanisms and bands. Such calculations were done in [82] by taking 
into account the large EPI interaction. By using the EPI spectral function a}{(o)F{(o) 
from tunnelling measurements and by assuming a}j.{(o)F{(o) a^{(o)F{(o) the authors 
have calculated A//((Oc, 2(0c) (and Az,(0, (Oc)) in the normal state and found a good agree- 
ment with experiments [79] - see Fig. 6. We stress that the recent elipsometric measure- 
ments of the dielectric function e{(o) [75] confirms this theoretical prediction [82]. 

From the above analysis we conclude that: (i) the interpretation of the partial sum- 
rule in HTSC oxides only in terms of the kinetic energy (or ©^^^(r)) is physically 
unjustified; (ii) the EPI interaction is strong and dominating scattering mechanism in 
the optical properties of the normal state. Reliable calculations of the partial sum-rule 
in the superconducting state are still missing, since in that case one should know much 
more details on the superconducting order A(k, co) and r;r(k, co), which are at present 
too ambitious task. 



A lot of experimental and theoretical works were devoted to the temperature de- 
pendence of resistivity p{T) in HTSC oxides. General properties of the resistivity in 
HTSC oxides are the following: (1) The resistivity is very anisotropic in single crys- 
tals where one has = {pc{T) / Pa-b{T)) » 1 at T above T^ - see [90], i.e. ^ 300 
in Lai.85Sro.i5CMC>4 and Nd\,'gsCeQ^isCuO^, rc ?a 20— 150 in Y Ba2Cu2,0'] -x, r^ ^ 10^ 
in Bi2Sr2CaCu20% depending also on the sample preparation, temperature etc. The 
anisotropy of the in-plane resistivity is much less, i.e. = {paa{T) / pt,b{T)) ~ 1 — 2, 
depending also on the sample preparation, temperature etc.; (2) The in-plane resistivity 
Pa-b{T) at room temperature is more than two orders of magnitude higher than that of 
the metallic Cm (where Pcu(Troom) ~ 1-5 JU^2cm),i.e. Pa-b(T) of HTSC oxides lies more 




(23) 



2.2.2. Resistivity p{T) 




FIGURE 6. Calculated theoretical T-dependence of high-energy part of the sum rule Ah{(Oc = 
1 .25eV, 2(0c = 2.5eV) by taking into account the electron-phonon interaction. The data from [82] 



in the semiconductor range and Pa-b{T) ^ PcuiJ); (3) Pa-biJ) r^T for T > Tc, which 
deviates at T > (800 — 1000) K and saturates at even higher temperatures, depending 
on samples etc.; (4) p^-^, varies from p^-biT) ~ T (with small residual resistivity) in 
optimally doped systems being pa-b{T) ~ T-^/^ in overdoped systems, as experiments 
on La2-xSrxCu04 show [89]; (5) In most samples of HTSC oxides the c-axis resistivity 
pc{T) shows a non-metallic behavior especially in samples with huge anisotropy along 
the c-axis, growing by decreasing temperature, i.e. {dpc{T) / dT) < 0, being supercon- 
ducting below 

We discuss briefly the in-plane resistivity pa-b{T) only, because its temperature 
behavior is a direct consequence of the quasi-2Z) motion of quasiparticles and of the 
inelastic scattering which they suffer. At present there is no consensus on the origin 
of the linear temperature dependence of the inplane resistivity Pa-b{T) in the normal 
state. As it is stressed several times many researchers are (erroneously) believing that 
such a behavior can not be due to the EPI? The inadequacy of this claim was already 
demonstrated by analyzing the dynamical conductivity a{(o). The inplane resistivity in 
HTSC oxides is usually analyzed by the Kubo approach, or by the Boltzmann equation. 
In the latter case p (r) is given by 

p{T)^^T,r{T) (24) 

7t CO 

TtriT) = - / d(0 . c4(6))F(6)), (25) 

i Jo sur((o/2T) 



where a}^{o))F {(o) is the EPI transport spectral function. It is well-known that at T > 
0£)/5 and for the Debay spectrum one has 

p{T)^^A^/^^^=p'T. (26) 

In HTSC oxides the reach and broad spectrum of a}^{(o)F {(o) is favorable for such 
a linear behavior. The measured transport coupling constant Xtr contains in principle 
all scattering mechanisms, although usually some of them dominate. For instance, the 
proponents of the spin-fluctuations mechanism assume that Xtr is entirely due to the 
scattering on spin fluctuations. However, by taking into account specificities of HTSC 
oxides the experimental results for the inplane resistivity pa-h{T) can be satisfactory 
explained by the EPI mechanism. From tunnelling experiments [37], [38], [39], [40], 
[41] one obtains that A 2 — 3 and if one assumes that Xtr ^ X and cOpi > (3—4) 
eV (the value obtained from the band- structure calculations) then Eq.(26) describes 
the experimental situation rather well. The plasma frequency cOpi which enters Eq.(26) 
can be extracted from optic measurements ((Opi^ex)^ i-e. from the width of the Drude 
peak at small frequencies. However, since ?itr ~ 0.25coj[{eV)p\iiacm/K) there is an 
experimental constraint on Xtr- The experiments [76] give that (Opi (2 — 2.5) eV and 
p' ^ 0.6 in oriented YBCO films, and p' 0.3 in single crystals of BISCO. These results 
makes a limit on Xtr ~ 0.9 — 0.4. 

So, in order to explain p{T) with small ?itr and high Tc (which needs large X) by 
the EPI it is necessary to have Xtr < {X/3). This means that in HTSC oxides the 
EPI is reduced in transport properties where ^ A. This reduction of cOp and Xtr 
means that they contain renormalization (with respect to the LDA results) due to various 
quasiparticle scattering processes and interactions, which do not enter in the LDA theory. 
In subsequent chapters we shall argue that the strong suppression of Xtr may have its 
origin in strong electronic correlations [17], [18], [19]. 

In conclusion, optic and resistivity measurements in normal state of HTSC oxides 
are much more in favor of the EPI than against it. However, some intriguing questions 
still remains to be answered: (i) which are the values of Xtr and cOpf. (ii) why one 
has Xtr ^ A: (iii) what is the role of the Coulomb scattering in a{(0) and p{T). 
The ARPES measurements (see discussion below) give evidence for the appreciable 
Coulomb scattering at higher frequencies, where r(tt)) pa Fq -|- AcCt) for co > co^ax with 
Ac ~ 0.4. So, in spite of the fact that the EPI is suppressed in transport properties it 
is sufficiently strong in order to dominate in some temperature regime. It may happen 
that at higher temperatures the Coulomb scattering dominates in p{T), which certainly 
does not disqualify the EPI as the pairing mechanism in HTSC oxides. For better 
understanding of p{T) we need a controllable theory for the Coulomb scattering in 
strongly correlated systems, which is at present lacking. 



2.3. Raman scattering in HTSC oxides 



If the elementary excitation involved in the Raman scattering are electronic we deal 
with the electronic Raman effect, while if an optical phonon is involved we deal with the 



phonon Raman effect. The Raman scattering in tlie normal and superconducting state 
of HTSC oxides is an important spectroscopic tool which gives additional information 
on quasiparticle properties - the electronic Raman scattering, as well as on phonons and 
their renormalization by electrons - the phonon Raman scattering. 



2.3.1. Electronic Raman scattering 

The Raman measurements on various HTSC oxides show a remarkable correlation 
between the Raman cross-section 5exp(<o) and the optical conductivity Oa-b{co), i.e. 

Sexpico) ~ [1 7,,(q) \^)FG)aa-b{G)), (27) 

where nfi((o) is the Bose function and 7ic(q) screened Raman vertex - see more in [2]. 
Previously it was demonstrated that Oa-tio)) depends on the transport scattering rate 
rtr{(0,T) where Ttrico^T) ~ T and ^^(co) T/co for (O < T, thus giving S(q, co) ^ 
Consti in that range. For CO > T one has coaa-bioi) ~ Const giving also 5(q, CO) ^ 
Const2. We have also demonstrated that the EPI with the very broad spectral function 
a^F{(o) (see Fig. 11 below) explains in a natural way (0,T dependence of <7a_^((o) 
and rfr(<y, r). So, the Raman spectra in HTSC oxides can be explained by the EPI in 
conjunction with strong correlations. This conclusion is supported by calculations of the 
Raman cross-section [91] which take into account the EPI with a}F{(o) extracted from 
the tunnelling measurements on YBa2Cu^0^j^x ^i^d Bi2Sr2CaCu20s+x [37]. They are in 
a good qualitative agreement with experimental results - see more in [2] . 

We stress again, that quite similar (to HTSC oxides) properties of the electronic 
Raman scattering (besides a((o), R{0)) and p{T)) were observed in experiments [92] 
on isotropic metallic oxides LaQ ^SrQ^CoO^ and CaQ ^Sro^RuO^ - see Fig. 7. To repeat 
again, in these compounds there are no signs of antiferromagnetic fluctuations (which 
are present in HTSC oxides) and the peculiar behavior is probably due to the EPI. 



2.3.2. Phonon Raman scattering 

Normal state. - The effect of the EPI on the Raman scattering is characterized by 
the Fano asymmetry parameter q{(0) - see more in [2]. If it is finite the line shape is 
asymmetric - the Fano effect (resonance). By decreasing q{o)) the phonon line shape 
becomes more asymmetric, which means stronger EPI (in case when continuum states 
are due to conduction carriers). The Fano resonance is experimentally found in HTSC 
oxide YBa2Cu30-]_s [94], where the line asymmetry is clearly seen in optimally doped 
(5 <g 1) systems, while it is absent in the insulating state (5 = 1). The existence of the 
Fano (asymmetric) line shape in HTSC oxides is a direct proof that the discrete phonon 
level interacts with continuum of states, which are conduction electrons in the metallic 
state - see Fig. 8. 
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FIGURE 7. Broad range Raman scattering spectra of Cao sSro ^RuOi, (broken line) and Lao ^Sro sCoO-i 
(solid line). Inset spectra of Tl2Ba2Ca2Cu30io, Bi2Sr2CaCu20%, YBa2Cu30j and Lai,%5Sro,t5Cu04. 
From [92]. 



Superconducting state. - It is well known that the renormalization of phonon fre- 
quencies and their life-times by superconductivity in LTSC materials is rather small - 
around one percent. The smallness of the effect is characterized by the parameter A/Ep 
which is very small in low temperature superconductors. However, A/Ef is much larger 
in HTSC oxides and already from that point of view one expects much stronger renor- 
malization effects. At the very beginning several Raman active phonon modes, with 
frequencies 128,153,333,437 and 501 cm~^, were detected in YBa2Cu20j and these 
modes are totally symmetric modes (with respect to the orthorhombic point group D2h)- 
(Icm^i = 29.9SGHz = 0.123985mey = 1.44A:) However, by using the approximate 
tetragonal symmetry (with the point group 1)4/,) the mode at (Osj^ = 333 cm^ trans- 
forms according to the Big representation, while the other modes according to the Ai^ 
one - see Fig. 9. The Fano resonance (asymmetric line shape) of the Big mode indicates 
an appreciable coupling of the lattice to the continuum, which in fact corresponds to 
the charge carries. It is interesting to note that the Aig modes in YBCO are weakly af- 
fected in the presence of superconductivity, while the Big mode softens by 9 cm~^ (by 
approximately 3 %) [93]. It is well established also that this softening is due to super- 




FIGURE 8. Fano resonance in YBa2Cu2,Oq-x. The asymmetry is seen for 1 12 and 337 cm phonons 
in the superconductor {x = 0). The semiconductor {x = 1) has Lorenzian line shapes. From [94]. 



conductivity and not due to, for instance, structural changes, because it disappears in 
magnetic fields higher than Hc2 . 

The frequency shift Scox and the phonon line width Ti in the superconducting state 
have been studied numerically in [97] for the case of the isotropic s — wave super- 
conducting gap (A(k) = A = const) and for strong coupling superconductivity. They 
have predicted the phonon-softening and line-width narrowing for ftJo < 2A, while for 
COo > 2A there is a phonon-hardening and line-width broadening. These predictions 
are surprisingly in agreement with experiments [93], in spite of the assumed isotropic 
s — wave pairing what is contrary to the experimentally well established d — wave pair- 
ing in YBCO. Later calculations of the renormalization of the fiig Raman phonon mode 
in the presence of the weak coupling d — wave superconductivity [100] show that if one 
assumes that cOb,^ < 2Aniax there is phonon softening accompanied with the line broad- 
ening below Tc. The latter is possible because of the gapless character (on a part of the 
Fermi surface) of J — wave pairing. In that respect the calculations of the phonon renor- 
malization based on the strong coupling d — wave superconductivity are of significant 
interest and still awaiting. 

Recent report [4] on the superconductivity-induced strong phonon renormaliza- 
tion of the Aig phonons at 240 and 390 cm~^ (by 6 and 18 % respectively) in 
HgBa2Ca^Cu^0 iQj^x {Tc = 123 K) - the so called Hg — 1234) compound, renders an 
additional evidence for the strong EPI in HTSC oxides. In [4] the EPI coupling constant 
is estimated to be rather large for the Aig phonons {^Aig ~ 0.08). Since there are 60 
phonon modes in HgBa2Ca^Cu/^0 \qj^x they are capable to produce large EPI coupling 
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FIGURE 9. Assignation of Ag modes according to calculations in [95]. In brackets are experimental 
phonon frequencies in cm~ ^ 1 15 cm~ ^ is the Ba mode, 150 cm~ ^ is the Cu2 mode,340 cm~ ^ {Big mode) 
and 445 cm~^ modes are due to vibration of 0(2, 3) ions in the Cu02 , while 505 cm~^ mode is due to 04 
ions. From [96]. 



constant A = Lv=i > 1 - see Fig. lO(a-b). A conservative estimation of the upper 
limit of Amax gives Amax ~ 60 X 0.08 = 4.8 wliich is, of course, far from the realistic 
value of A < 2. In any case this analysis confirm that the EPI of some Raman modes 
in HTSC oxides is strong. To this point, very recent Raman scattering measurements 
on the (Ca,C) — 1234 compound with 7^ = 117 ^ reveal strong superconductivity 
induced phonon self-energy effects [5]. The Ai^ phonons at 235 cm~^ and 360 cm~^ 
(note cOpi < 2Ao), which involve vibrations of the plane oxygen with some admixture 
of Ca displacements, exhibit pronounced Fano line shape (in the normal and supercon- 
ducting state) with the following interesting properties in the superconducting state: 
(i) the phonon intensity is increased substantially; (ii) both phonons soften; (Hi) the 
phonon line width (of both phonons) increases dramatically below T^ passing through a 
maximum slightly below T^, and decreases again at low T but remaining broader than 
immediately below Tc. This line broadening is difficult to explain by s — wave pairing, 
where the line narrowing is expected, but it can be explained by superconducting pairing 
with nodes in the quasiparticle spectrum, for instance by J — wave pairing [98]. The 
large EPI coupling constants for these two modes are estimated from the asymmetric 
Fano line shape, which gives A235 = 0.05 and A360 = 0.07 (note in YBCO Aaj^ = 0.01 
for (OAig = 440 cm~^ and Xs^g = 0.02 for (Ogj^ = 340 cm~^, rather small values) giving 
the upper value for the total coupling constant Amax = 4. This result gives additional 
important evidence for the strong EPI in HTSC oxides. 



2.3.3. Electron-phonon coupling in Raman scattering 

We would like to stress the importance of the (phonon) Raman scattering measure- 
ments for the theory of the EPI in HTSC oxides. The covalent part of the EPI is due to 
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FIGURE 10. The fitted frequency cOp, line-width Tp, asymmetry parameter q, and the phonon intensity 
Ip of the Hg — 1234 Raman spectra in the Aig mode measured in xW polarization with 647.1 nm laser 
line: (a) at 240 cm~^; [b) at 390 cm~^ - from [4]. 



the Strong covalency of the Cu and O orbitals in the Cu02 planes. In that case the EPI 
coupling constant is characterized by the parameter ("field") E'^'^ ~ dtp-d /dR ~ qotp-d, 
where tp-^ is the hopping integral between Cu{d^2_y2) and 0{pxj) orbitals and the 

length q^^ characterizes the spacial exponential fall-off of the hopping integral tp^^- 
The covalent EPI is unable to explain the strong phonon renormalization (the self-energy 
features) in the Big mode in YBazCu^Oj by superconductivity, since in this mode the 
O-ions move along the c — axis in opposite directions, while for this mode dtp-d/dR is 
zero in the first order in the phonon displacement. Therefore the EPI in this mode must 
be due to the ionic contribution to the EP interaction which comes from the change 
in the Madelung energy as it was first proposed in [47], [48]. Namely, the Madelung 
interaction creates an electric field perpendicular to the Cu02 planes, which is due to 
the surrounding ions which form an asymmetric environment. In that case the site en- 
ergies ef contain the matrix element e^"" = (v^/ | 0(r) | where | i/Zj) is the atomic 
wave function at the /-th site, while the potential ^ (r) steams from surrounding ions. In 
simple and transition metals the surrounding ions are well screened and therefore the 
change of e'"" in the presence of phonons is negligible, contrary to HTSC oxides which 
are almost ionic compounds (along the c-axis) where the change of e'"" is appreciable 



and characterized by the field strength E""^ = V /un- Here, V is the characteristic poten- 
tial due to surrounding ions and a„ is the distance of the neighboring ions. Immediately 
after the discovery of HTSC oxides in many papers [1 14], [60], [99] it was (incorrectly) 
assumed that the covalent part dominates the EPI in these materials. The calculation of 
Tc by considering only covalent effects [114], [60] gave rather small (~ 10 — 20 ^ 
in YBCO, and 20 — 30 in Lai_%sSrQ_\5CuO^). It turns out that in HTSC oxides the op- 
posite inequality E'^" ^ E^""^ is realized for most c-axis phonon modes, on which basis 
the renormalization of the Raman B\g mode can be explained - see more in [2]. This is 
supported by detailed theoretical studies in for the YBCO compound [48], [49], where 
it is calculated the change in the ionic Madelung energy due to the out of plane oxygen 
vibration in the B\g mode. Similarly as in YBCO, the large superconductivity-induced 
phonon self-energy effects in H gBa2CaT,Cu^0 iq^^ and in {Cu,C)Ba2Ca'iCu^0\Q+x for 
the Aig modes are also due to the ionic (Madelung) coupling. In these modes oxygen 
ions move also along the c — axis and the ionicity of the structure is involved in the 
EPI. This type of the {long-range) EPI is absent in usual isotropic metals (LTSC super- 
conductors), where the large Coulomb screening makes it to be local. Similar ideas are 
recently incorporated into the Eliashberg equations in [133], [132]. The weak screening 
along the c-axis, which is due to the very small hopping integral for carrier motion, is 

reflected in the very small plasma frequency £0^'^'' along this axis. Since for some optical 

phonon modes one has (Oph > (Op ^ then nonadiabatic effects in the screening are impor- 
tant. The latter can give rise to much larger EPI coupling constant for this modes [52], 
[53]. 

In conclusion, the electron and phonon Raman scattering measurements in the normal 
and superconducting state of HTSC oxides give the following important results: (a) 
phonons interact strongly with the electronic continuum, i.e. the EPI is substantial; (b) 
the ionic contribution (the Madelung energy) to the EPI interaction for c-axis phonon 
modes gives substantional contribution to the (large) EPI coupling constant (A > 1). 



2.4. Tunnelling spectroscopy in HTSC oxides 

Tunnelling methods are important tools in studying the electronic density of states 
N{q)) in superconductors and in the past they have played very important role in investi- 
gating of low Tc-superconductors. By measuring the current- voltage (/ — V) characteris- 
tic in typical tunnelling junctions (with large tunnelling barrier) it was possible from the 
tunnelling conductance G{V){= dl/dV) to determine N{o}) and the superconducting 
gap as a function of temperature, magnetic field etc. Moreover, by measuring of GiV) at 
voltages > A in the NIS (normal metal - isolator - superconductor) junctions it was 
possible to determine the Eliashberg spectral function a^F{co) (which is due to some 
bosonic mechanism of quasiparticle scattering) and finally to confirm (definitely) the 
phonon mechanism of pairing in LTSC materials, except maybe heavy fermions [101]. 
We shall discuss here only the results for a^F{co) obtained from I — V measurements, 
while a more extensive discussion of other aspects is given in [2]. 



2.4.1. I — V characteristic and a^F{(o) 



If one considers a NIS contact where the left (L) and right (R) banks of the contact can 
be normal (A'^) metal or superconductor (5), respectively, with very small transparency 
then tunnelling effects are studied in the framework of the tunnelling Hamiltonian 
Ht = Lk,p(^k,pC|^^Cq^ + h.c). In that case the single-particle tunnelling current is given 
by the formula [102] 

/,p(y) = 2.£|rk,p|2x 

k,p 

/OO 
d(0AN{k, (o)As{-p,(0 + eV)[nf{(o)-nf{(0 + eV)]. (28) 
-CXD 

The single-particle spectral function AjY(5)(k, to) is related to the imaginary part of 
the retarded single particle Green's function, i.e. A(k, (o) = —ImG^^^ilL, (£>)/%, while 
the tunnelling matrix element | Tk p p is derived in the quantum-mechanical theory 
of tunnelling through the barrier - see [2]. Note, in the superconducting state A(k, ca) 
depends on the superconducting gap function A(k, (o), which is on the other hand a 
functional of the spectral function oc^F(tt)). The fine structure in the second derivative 
cP-I/dV^ at voltages above the superconducting gap is related to the spectral function 
a^F{co). For instance, plenty of break-junctions made from Bi — 2212 single crystals 
[37] show that negative peaks in d^I/dV^, although broadened, coincide with the peaks 
in the generalized phonon density of states Gph{0)) measured by neutron scattering - see 
more in [2]. Note, the reported broadening of these peaks might be partly due tod— wave 
pairing in HTSC oxides. The tunnelling density of states Nt{V) ~ dl/dV shows a gap 
structure and it was found that 2A/Tc = 6.2 — 6.5, where Tc = 14 — S5K and A is some 
average value of the gap. By assuming s — wave superconductivity [37] and by solving 
the MR problem (inversion of Eliashberg equations), the spectral function a^F{(x)) is 
obtained which gives A ^ 2.3. Note, in extracting A [37] the standard value of the 
effective Coulomb parameter /i* ^ 0. 1 is assumed. Although this analysis [37] was done 
by assuming s — wave pairing it is qualitatively valuable procedure also in the case of 
d — wave pairing, because one expects that d — wave pairing does not spoil significantly 
the global structure of d^I/dV^ at eV > A, but introducing mainly a broadening of peaks. 
The latter effect can be partly due to an inhomogeneity of the gap. The results obtained 
in [37] were reproducible on more than 30 junctions, while in 5/(2212) — GaAs and 
5/(2212) —Am planar tunnelling junctions similar results were. Several groups [39], [40], 
[41] have obtained similar results for the shape of the spectral function a^F{(o) from 
the I — V measurements on various HTSC oxides as shown in Fig. 11. The results shown 
in Fig. 1 1 leave no much doubts on the effectiveness of the EPI in pairing mechanism 
of HTSC oxides. In that respect recent tunnelling measurements on Bi2Sr2CaCu20^ 
[103] are impressive, since the Eliashberg spectral function a^F{(o) was extracted from 
the measurements of (P-ljdY^ . The obtained a^F(£o) has several peaks in the broad 
frequency region up to 80 meV - see Fig. 1 1 (curve Shimada et al.), which coincide rather 
well with the peaks in the phonon density of states F(ft)). Moreover, the authors of [103] 
were able to extract the coupling constant for modes laying in (and around) these peaks 
and their contribution to T^ . They managed to extract the EPI coupling constant, which 
is unexpectedly very large, i.e. 'k{= 2 / d(£)a}F{(i>)J (o) = fti 3.5. Since almost all 
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FIGURE 11. The spectral function a?'F{(o) obtained from measurements of G(y) by various groups 
on various junctions: Vedeneev et al. [37], Gonnelli et al. [41], Miyakawa et al. [39], Shimada et al.[38]. 
The generalized density of states GPDS for Billll is plotted at the bottom - from [38]. 



phonon mode contributes to A, this means that on the average each particular phonon 
mode is moderately coupled to electrons thus keeping the lattice stable. Additionally, 
they have found that some low-frequency phonon modes corresponding to Ca, Sr and 
Ca vibrations are rather strongly coupled to electrons, similarly as the high frequency 
oxygen vibrations along the c-axis do. These results confirm the importance of the axial 
modes in which the change of the Madelung energy is involved, thus supporting the idea 
conveyed through this article of the importance of the ionic Madelung energy in the EPI 
interaction of HTSC oxides. 

In conclusion, the common results for all reliable tunnelling measurements in HTSC 
oxides, including Bai-xKxBiOT, too [104], [16], is that no particular mode can be singled 
out in the spectral function a^F{(o) as being the only one which dominates in pairing 
mechanism. This important result means that the high T^ is not attributable to a particular 
phonon mode in the EPI mechanism, since all phonon modes contribute to X. Having in 
mind that the phonon spectrum in HTSC oxides is very broad (up to 80 meV ), then the 
large EPI coupling constant (A 2) in HTSC oxides is not surprising at all. We stress. 



that compared to neutron scattering experiments the tunnelling experiments are superior 
in determining the EPI spectral function a^F((o). 



2.5. Isotope effect in HTSC oxides 

The isotope effect has played an important role in elucidating the pairing mechanism 
in LTSC materials. Note, the standard BCS theory predicts that for the pure phonon- 
mediated mechanism of pairing the isotope coefficient a = —dXwTc/dlviM, where M is 
the ionic mass, takes its canonical value a = 1/2. However, later on it was clear that 
a can take values less (even negative) then its canonical value in the phonon-mediated 
mechanism of pairing if there is pronounced Coulomb pseudopotential jU* - see more in 
[2]. 



2.5.1. Experiments on the isotope coefficient a 

A lot of measurements of Uq and Ocu were performed on various hole-doped and 
electron-doped HTSC oxides and we give a brief summary of the main results [105]: 
(1) The O isotope coefficient ao strongly depends on the hole concentration in the hole- 
doped materials where in each group of HTSC oxides (YBa^Cu^tO-j-x, or La2-xSryCuO^ 
etc.) a small oxygen isotope effect is observed in the optimally doped (maximal Tc) sam- 
ples. For instance ao ~ 0.02 - 0.05 in YBa2Cu30j with Tc^max 91 K, ao ^0.1-0.2 in 
Lai 855ro.i5CM04 with Tc^umx ~ 35 K; ao ~ 0.03 — 0.05 in Bi2Sr2CaCu20^ with T^max ~ 
76 K; ao ~ 0.03 and even negative (—0.013) in Bi2Sr2Ca2Cu20\Q with max ~ llO/sT; 
the experiments on Tl2Can_iBaCun02n+A = 2,3) with Tcmax ~ 121 K are still un- 
reliable and ao is unknown; ao < 0.05 in the electron-doped {Ndi-xCex)2Cu04 with 
7c,max ~ 24 K. (2) For hole concentrations away from the optimal one, Tc decreases 
while ao increases and in some cases reaches large value ao ~ 0.5 - see Fig. 12 for La 
compounds. This holds not only for parent compounds but also for systems with substi- 
tutions, like {Yi_x-y^^xCay)Ba2Cu30j, Yi_yCayBa2Cu404 and Bi2Sr2Ca\-xYxCu20^. 
Note, the decrease of T^. is not a prerequisite for the increase of ao- This became clear 
from the Cu substituted experiments YBa2{Cui-xZnx)30j where the decrease of T,. (by 
increasing of the Zn concentration) is followed by only small increase of ao [106]. Only 
in the case of very low <20 K then ao becomes large, i.e. ao> OA. (3) The largest 
ao is obtained even in the optimally doped compounds like in systems with substitu- 
tion, such as Lai.85S'ro.i5CMi_xMx04, M = Fe, Co, where ao^l-3 for x ^ 0.4 %. (4) In 
La2-xMxCu04 there is a Cu isotope effect which is of the order of the oxygen one, i.e. 
acu ~ OCo giving acu + cuo ~ 0.25 — 0.35 for optimally doped systems (x = 0. 15). In the 
case when x — 0. 125 with Tc -C r^^max one has acu ~ 0.8 — 1 with acu + ao^ 1.8. The 
appreciate copper isotope effect in La2-xMxCu04 tells us that vibrations of other than 
oxygen ions could be important in giving high T^. The latter property is more obvious 
from tunnelling measurements, which are discussed above. (5) There is negative Cu iso- 
tope effect in the oxygen-deficient system YBa2Cu30j-x where acu is between —0.14 
and —0.34 if T^ lies in the 60 K plateau. (6) There are reports on small negative ao in 
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FIGURE 12. The oxygen isotope exponent Uq for: (a) La2~xSrxCu04 as a function of Sr con- 
centration - from [105]. The oxygen isotope exponent ao as a function of T^ for: (b) YBa2Cu^0-i. 
1: {Yx_x^Vx)Ba2CuiOr, 2: YBa2-xLa^Cu-iOr, 3: YBa2{Cui-xCox)iOr, 4: Y Ba2{Cui-x2nxhOr, 5: 
YBa2{Cui-xFex)^0-i. (c) Lai,^sSrQ,isCu04. 1: Lai.g55ro.i5(CMi_;cM;c)04; 2: iai.855ro.i5(CMi_;cZn;t)04; 
3: Lai.85-S'ro.i5(CMi_;cCo;c)04; 4: Lai,%5SrQ,i5{Cui-xFex)OA - from [105]. 



some systems like YSr2CuzO'] with ao -0.02 and in B/SCO - 2223 {Tc = 110 ^) 
where ao ~ —0.013 etc. However, the systems with negative ao present considerable 
experimental difficulties, as it is pointed out in [105]. 

The above enumerated results, despite experimental difficulties, are more in favor 
than against of the hypothesis that the EPI interaction is strongly involved in the pairing 
mechanism of HTSC oxides. By assuming that the experimental results on the isotope 
effect reflect an intrinsic property of HTSC oxides one can rise a question: which 
theory can explain these results? Since at present there is no consensus on the pairing 
mechanism in HTSC materials there is also no definite theory for the isotope effect. 
Besides the calculation of the coupling constant X any microscopic theory of pairing 
is confronted also with the following questions: (a) why is the isotope effect small in 
optimally doped systems and (b) why a increases rapidly by further under(over)doping 
of the system? 

It should be stressed, that at present all theoretical approaches are semi-microscopic, 
but what is interesting most of them indicate that in order to explain the rather unusual 
isotope effect in HTSC materials one should invoke iht forward scattering peak in the 
EPI [2]. 

In conclusion, experimental investigations of the isotope effect in HTSC oxides have 
shown the importance of the EPI interaction in the pairing mechanism. 



2.6. ARPES experiments in HTSC oxides 



2.6.1. Spectral function A , co ) from ARPES 

The angle-resolved photoemission spectroscopy (ARPES) is nowadays a leading 
spectroscopy method in the solid state physics. The method consists in shining light 
(photons) with energies between 20 — 1000 eV on the sample and by detecting momen- 
tum (k)- and energy ((j!))-distribution of the outgoing electrons. The resolution of ARPES 
is drastically increased in the last decade with the energy resolution of f« 2 meV (for 
photon energies ~ 20 eV) and angular resolution of A0 ^ 0.2°. The ARPES method is 
surface sensitive technique, since the average escape depth {Igsc) of the outgoing elec- 
trons is of the order of lesc ~ 10 A. Therefore, one needs very good surfaces in order that 
the results be representative for the bulk sample. In that respect the most reliable stud- 
ies were done on the bilayer Bi2Sr2CaCu20% {Billll) and its single layer counterpart 
Bi2Sr2CuO() (5/2201), since these materials contains weakly coupled BiO planes with 
the longest interplane separation in the HTSC oxides. This results in a natural cleavage 
plane making these materials superior to others in ARPES experiments. After a dras- 
tic improvement of sample quality in others families of HTSC materials, became the 
ARPES technique a central method in theoretical considerations. Potentially, it gives in- 
formation on the quasiparticle Green's function, i.e. on the quasiparticle spectrum and 
life-time effects. The ARPES can indirectly give information on the momentum and en- 
ergy dependence of the pairing potential. Furthermore, the electronic spectrum of the 
HTSC oxides is highly quasi-2D which allows an unambiguous determination of the 
momentum of the initial state from the measured final state momentum, since the com- 
ponent parallel to the surface is conserved in photoemission. In this case the ARPES 
probes (under some favorable conditions) directly the single particle spectral function 
A(k, (o). 

In the following we discuss only those ARPES experiments which give us evidence 
for the importance of the EPI in HTSC oxides - see detailed reviews in [107], [108]. 

The photoemission measures a nonlinear response function of the electron system, 
since the photo-electron current (j(l)) at the detector is proportional to the incident 
photon flux (square of the vector potential A), i. e. schematically one has 

(j(l)) ^ (j(2)j(l)j(3))A(2)A(3), (29) 

and integration over bar (1 = (x,0 indices is understood. The correlation function 
(j(2)j(l)j(3)) describes all processes related to electrons, such as photon absorption, 
electron removal and electron detection, are treated as a single coherent process. In this 
case the bulk, surface and evanescent states, as well as surface resonances should be 
taken into account - the so called one-step model. 

Under some conditions the one-step model can be simplified by an approximative, but 
physically plausible, three-step model. In this model the photoemission intensity 

Itot{KO))^l-l2-h (30) 

is the product of three independent terms: (1) / - describes optical excitation of the 
electron in the bulk; (2) I2 - the scattering probability of the travelling electrons; (2) 73 - 



the transmission probability through the surface potential barrier. The central quantity in 
the three-step model is /(k, co). To calculate it one assumes the sudden approximation, 
i.e. that the outgoing electron is moving so fast that it has no time to interact with the 
photo-hole - see more in [107], [108]. It turns out that 7(k, ft)) can be written in the form 
[107], [108] (fork = k||) 

/(k, CO) ~ /o(k, u)/(ft))A(k, ft)). (31) 

/o(k, v) ~| ( V^/ I pA 1 1//,) p where | pA 1 1//,) is the dipole matrix element and depends 
on k, polarization and energy v of the incoming photons, /(ft)) = 1/(1 +exp{Q)/T}) is 
the Fermi function and A(k, ft)) = — ImG(k, ft))/;r is the quasiparticle spectral function. 
In reality because of finite resolution of experiments, in k and ft), /(k, ft)) should be 
convoluted by the ft)-convolution function R{(o) and k-convolution function 2(k). It 
must be also added the extrinsic background B, which is due to secondary electrons 
(those which escape from the sample after having suffered inelastic scattering events 
coming out with reduced kinetic energy). 

By measuring A(k, ft)) one can determine L(k, ft)) = Li (k, ft)) -|- i'L2{K ®) 

MK = _ ^^(k) - Li (k, ft))]2 + p2(k, ft))]2 • ^^^^ 

^o(k) = Ck — jU is the bare quasiparticle energy. For instance in the case of the Landau- 
Fermi liquid A (k, ft)) can be separated into the coherent and incoherent part 

A(k, ft)) = Zk _ ^ ^^^^^ ^ +Ainch{K CO), (33) 

where = 1/(1 — dLi/dco), t,{k) = Zk((^o(k) +£i) and Fk = Zk | E2 | calculated at 
ft) = (k). For small ft) one has ^ (k) »| ^2 | and Fk [(TtT)^ + <^^(k)]. 

In some period of the HTSC era there were a number of controversial ARPES results 
and interpretations, due to bad samples and to the euphoria with exotic theories. For 
instance, a number of (now well) established results were questioned in the first ARPES 
measurements, such as: the shape of the Fermi surface, which is correctly predicted by 
the LDA band-structure calculations; bilayer splitting in B/2212, etc. 

We summarize here important ARPES results which were obtained recently, first in 
the normal state [107], [108]: (Nl) There is well defined Fermi surface in the metallic 
state - with the topology predicted by the LDA; (N2) the spectral line are broad with 
I Z2(k, ft)) I ft) (or ~ r for r > ft)); (N3) there is a bilayer band splitting in B/2212 (at 
least in the overdoped state); (N4) at temperatures Tc < T < T* and in the underdoped 
HTSC oxides there is a d-wave like pseudogap Apg(k) ~ Apgfi{coskx — cosky) in the 
quasiparticle spectrum; (N5) the pseudogap A^^ increases by lowering doping; (N6) 
there is evidence for the strong EPI interaction and characteristic phonon energy (Oph - 
dXT >Tc. Namely, in all HTSC oxides which are superconducting there are kinks in the 
quasiparticle dispersion in the nodal direction (along the (0,0) — {n,n) line) at around 

aP^^ ~ (60 - 70) meV [43] - see Fig. 13, and around the anti-nodal point (;r,0) at 40 
meV [109] -see Fig. 14. 




FIGURE 13. Quasiparticle dispersion of Bi22l2, B/2201 and LSCO along the nodal direction, plotted 
vs the momentum k for (a) — (c) different doings, and (d) — (e) different T; black arrows indicate the kink 
energy; the red arrow indicates the energy of the q ~ (7^,0) oxygen stretching phonon mode; inset of (e)- 
T-dependent E' for optimally doped Bi2212; (/) - doping dependence of A' along (0,0) — {n, n) for the 
different HTSC oxides. From [43] 



In the superconducting state ARPES results are the following [107], [108]: (SI) there 
is an anisotropic superconducting gap in most HTSC compounds, predominately of d- 
wave like, A5c(k) ~ Ao(cosfc,: — cosfcy) with /Tc ~ 5 — 6; (S2) the dramatic changes 
in the spectral shapes near the point (7r,0), i.e. a sharp quasiparticle peak develops at 
the lowest binding energy followed by a dip and a broader hump, giving rise to the so 
called peak-dip-hump structure; (S3) the kink at (60 — 70) meV is surprisingly unshifted 
in the superconducting state -[43]. To remind the reader the standard Eliashberg theory 

the kink should be shifted to cOph + Aq. (S4) the anti-nodal kink at (0^^^^ ~ 40 meV is 

shifted in the superconducting state by Aq, i.e. (o'j^^^ ^^ph^ +^o = (65 — 70)meV since 
5 = (25-30)meV-see[109]. 




FIGURE 14. Quasiparticle dispersion E{k) in the normal state (al, bl, c), at 107 K and 115 K, along 
various directions (p around the anti-nodal point. The kink at £ = AOmeV is shown by the horizontal arrow. 
(a2 and b2) is E{k) in the superconducting state at 10 K with the shifted kink to IQmeV. (d) kink positions 
as a function of (/) in the anti-nodal region. From [109] 

2.6.2. Theory of the ARPES kink 

We would like to point out that the breakthrough-experiments done by the Shen group 
[43], [109] shown in Fig. 13, Fig. 14, with the properties (N6), (S3) and (S4) - which 
we call the ARPES shift-puzzle, are the smoking-gun experiments for the microscopic 
theory of HTSC oxides. Namely, any theory which reflects to explain the pairing in 
HTSC oxides must solve the shift-puzzle. 

In that respect the recent theory [44], which is based on the existence of the forward 
scattering peak (which is due to strong correlation in the EPI) in the EPI - the FSP 
model, was able to explain this puzzle in a consequent way. The FSP model (see more 
in [2], [44]) contains the following basic ingredients: (i) the electron-phonon interaction 
is dominant in HTSC and its spectral function a^F(k,k', D.) ^ a^F{(p, (p' , ^) {(p is the 
angle on the Fermi surface) has a pronounced forward scattering peak due to strong 
correlations. Its width is very narrow | k — k' 1^-^ kf even for overdoped systems [17], 
[18], [19]. In the leading order a^F{q),(p' ~ d{(p — (p'); (ii) the dynamical part 
(beyond the Hartree-Fock) of the Coulomb interaction is characterized by the spectral 
function 5c (k, k', Q). The ARPES shift puzzle implies that 5c is either peaked at small 
transfer momenta | k — k' | , or «Y is so small that the shift is weakly affected and is beyond 
the experimental resolution of ARPES. We assume that the former case is realized; (Hi) 



The scattering potential on non-magnetic impurities has pronounced forward scattering 
peak, which is also due to strong correlations [17], [18], [19]. The latter is characterized 
by two rates Yi(2)- The case /i = 72 mimics the extreme forward scattering, which does 
not affect pairing. On the other hand , 72 = describes the isotropic exchange scattering 
- see discussion in 

The Green's function is given hy G^— 1/ {i(Ok — — 2^yk(tt))) = — + <^k) / (w^ + 
+ A|) where in the k = (k, ft)). In the FSP model the equations for 0% and A^^. are [44] 

^ X-i m{n — m)6)mm r 
ft^,^ = ft^ + ;rr £ ^ "''^ + E^,^ , (34) 

K,<p = TiTl^ 'P ^ +Agy, (35) 

where 

'^l(2),(p(« -m)= Xph,(p{n -m) + 4jn7i(2),<p 
with the electron-phonon coupling function 



f°° i2 
Xph,cpin) = 2 daal,^^F^{a)^^-^. (36) 



Since the EPI and in Eq.(34-35) has a local form as a function of the angle (p, 
then the equation for cbn^q, has also local form, which means that the different points 
on the Fermi surface are decoupled. In that case G)n,q) depends on the local value of 

the gap A„ (p pa AqCOsI^. Just this property is important in solving the ARPES shift 
puzzle. So, in the nodal point (cp = 71 /4) one has A„ .,p = and the quasiparticle spectrum 
given by £■ — — Eyt(^:^M,(p = 0) = is unaffected by superconductivity, i.e. the kink 
is unshifted. This is exactly what is seen in the experiment of the Shen group [43] - see 
Fig. 13. In the case of the antinodal point (<p ~ 7t/2) there is a singularity at AOmeV'm the 
quasiparticle spectrum (i^sing) in the normal state - see Fig. 14. The analytic and numeric 
calculations of Eq.(34) show that this singularity is shifted by Aq in superconducting 
state, i.e. £'sing £^sing +-^0- This is exactly what is seen in the recent experiment on 
BISCO [109] - see Fig. 14, where the singularity of the normal state spectrum at 40 meV 
is shifted to (65 — 70) meV in the superconducting state, since Aq ~ (25 — 30) meV . The 
FSP model explains in the natural way also the peak-dip-hump structure in A(k, ft)) - for 
more details see [44]. 



2.6.3. ARPES and the EPI coupling constant X 

One can rise the question - is it possible to extract the coupling constant X from 
ARPES measurements. As we have seen above, by assuming that the three-step model 
holds, where /(k, ft)) ~ A(k, ft)), then one possibility is by measuring the kink in the 
quasiparticle renormalization, i.e. by measuring the real part of the self-energy Ei (k, ft)). 



These measurements [43], [109] give A^);^^^ ~ 1 in both nodal and antinodal direction. 
Another possibility is by measuring the width (Akfwico)) of the momentum distribution 
curves (MDCs) which give the imaginary part E2(ft), T) via 

i:2ico,T)^^VFAkFwico). (37) 

In that respect very indicative are recent measurements [45] of E2((U,r) around 
the nodal point in a number of HTSC compounds, such as the superstructure 
free Bh-xPbxSr2CaCu20^^s (Bi{Pb) - 2212), Bi2Sr2CaCu20s+s (Bi - 2212) and 
Bi2Sr2^xLaxCu20^j^S (Bi — 2201). In the analysis of ARPES spectra the authors in 
[45] have assumed that there are two CQ-dependent contributions to E - the Fermi liquid 
contribution Lfi and the part due to the interaction via bosonic excitations (let say 
phonons and spin-fluctuations), i.e. E2 = E2,fl + E2,fi + E2,;mp - see Fig. 15. We stress, 
that in [45] it is assumed that vp = AeV A. 

In the framework of this procedure the theoretical analysis [46] gives small coupling 

^ARPES ^ which is extracted from the slope of E2 in the interval 0.05 eV < (O <0.\ 
eV - see the gray solid line in Fig. 15b). Such a small coupling gives very small and 

(2) 

none of pairing mechanisms is effective. Furthermore, the small value of ^\^p^g (most 

(2) 

measurements give ^^^^rpes < ^•^) " which is extracted from E2, is a generic property 
of most ARPES measurements. This means, that we are confronted with a trilemma: 
(1) to abandon the boson-fermion separation procedure of E((o) done in [45], (2) to 
abandon the Eliashberg theory, (3) to abandon the interpretation of ARPES data within 

the three-step model. 

In some sense the situation in ARPES measurements with A-^^p^^ <C X resembles 
the one in transport measurements where Aj^ <^ A, i.e. . This problem deserves further 
investigation. 



3. EPI IN HTSC OXIDES 

In the following we present briefly some elements of the general theory of the strong 
EPI and its low-energy version. In the latter, high-energy processes are integrated out 
and the low-energy phenomena are governed by the high-energy vertex functions r^, the 
excitation potential Eq (part of the self-energy due to the Coulomb interaction) and the 
EPI coupling constants gEP,ren ~ - see more in [2]. However, this procedure was never 
performed in its full extent, because of difficulties to calculate Eq and . Therefore, the 
EPI coupling constant was as a rule calculated by some other methods. Usually in LTSC 
materials the EPI is calculated by using the local-density functional (LZM) method, 
which is suitable for ground state properties of crystals (matter) and which is based on 
an effective electronic crystal potential Vg. Since in principle Vg may significantly differ 

from Eq then the LDA calculated coupling constant g^pp'^^ can be also very different from 

the real coupling constant gEp. The calculation of g^^p^^ complicated, even in the LDA 
method, and further approximations are necessary, like for instance the rigid-ion (RI) and 




Energy (eV) 

FIGURE 15. T- and o-dependence of E2 for the nodal quasiparticles in optimally doped Bi(Pb)-2212. 
(a) - the full width at half maximum of the ARPES intensity. The gray solid line is the Fermi liquid 
parabola obtained by fitting the data for highly overdoped sample (OD69) at 130 K (see inset), (b) - the 
bosonic part E2.B for various T.From [45] 



rigid muffin-tin (RMTA) approximations. These approximation were justified in simple 
metals. However, these approximations are inadequate for HTSC oxides, because they 
fail to take into account correctly the long-range forces (Madelung energy - see below) 
and strong electronic correlations. Strictly speaking the EPI does not have meaning in 
the LDA method - see more in [2], because the latter treats the ground state properties of 
materials, while the EPI is due to excited states and inelastic processes in the system. 



3.1. General strong coupling theory of the EPI 



It is based on the fully microscopic electron-ion Hamiltonian for the interacting elec- 
trons and ions in a crystal - see for instance [110], [111], [11 2], and comprises electrons 
interacting between themselves as well as with ions and ionic vibrations. In order to 
describe superconductivity the Nambu-spinor V^^(r) is introduced which operates in the 
electron-hole space ^^'^(r) = {y/^{r) Yli''')) (analogously for the column ^^(r)) with 
VA|(r), V>| (r) as annihilation and creation operators for spin up, respectively etc. The 
microscopic Hamiltonian which in principle should describe the normal and supercon- 
ducting state of the system contains three parts H = Hg + Hi + Hg-i- The electronic 
Hamiltonian Hg, which describes the kinetic energy and the Coulomb interactions of 
electrons, is given in the second-quantization by 

He = I d^rxi/\r)T3eo{p)xi/{r)+ 

+i I jVjV>t(r)f3v/(r)K(r-r')v/^(r')f3V^(r'), (38) 

where eo{p) — p^/2m is the kinetic energy of electron and Vc{r — r') — e^/ \ r — r' \ 
is the electron-electron Coulomb interaction. Note, that in the electron-hole space the 
pseudo-spin (Nambu) matrices f,-, / = 0, 1, 2, 3 are Pauli matrices. 

The lattice Hamiltonian (describes lattice vibrations Uan of ions enumerated by n) is 
given by 

^ n ^ n,m,a ^ n,m 

+ 1 ("a«-"am)(%«-%JV«V^V,(RO-RO)+Hr^ (39) 

n,m,a,p 

The first term in Eq.(39) is the kinetic energy of vibrating ions (with charge Ze), 
ViiCR^ — R|]J = Z^e^ j I R|J— RJ^ I is the bare ion-ion interaction in equilibrium, while 
the third and fourth terms describe the change of Ya by lattice vibrations with the ion- 
displacement is = R,j — R|]. The term describes higher anharmonic terms with 

respect to . The theory which we describe below holds for any kind of anharmonicity. 

The electron-ion Hamiltonian describes the interaction of electrons with the equilib- 
rium lattice and with its vibrations, respectively 

He-i = Y, I d\Ve-i{r-Kl)^f\r)T3^f{r)+ I d\^{r)^f\r)T3^f{r), (40) 
4>(r) = £[K-,-(r -K- fin) - ye-iir - R°) . 

n 

Here, V'e_,(r — R^) is the electron-ion potential - see [2]. The second term which de- 
pends on the lattice distortion operator <I>(r) describes the interaction of electrons with 
harmonic (~ ««„) (or anharmonic ~ m^„, ^: = 2, 3...) lattice vibrations. 



Based on the above Hamiltonian one can in principle calculate the electron and 
phonon Green's functions 

G(l,2) = -(rv/(l)v/t(2)) (41) 

and 

D(l-2) = -(r4>(l)4>(2)), (42) 
respectively. The solution of these equations is written in the form of Dyson's equations 

(5-1(1,2) = Go 1(1,2) -i:(l, 2) (43) 

and 

5-1(1,2) =D-i(l,2)-n(l, 2), (44) 

where Go 1(1,2) and 1(1,2) are the bare inverse electron and phonon Green's func- 
tion, respectively. The nontrivial effects of interactions are hidden in the self -energies 
L(l,2) and n(l,2). Here, 1 = (ri,Ti), where Ti is the imaginary time. The calculation 
of E is simplified by using the Migdal adiabatic approximation [139], which incorpo- 
rates the experimental fact that in most metals the characteristic phonon (Debye) energy 
of lattice vibrations (Ou is much smaller than the characteristic electronic Fermi energy 
Ep {(Od <^ Ep). Using this fact Migdal formulated a theorem which claims that in the 
self-energy L one should keep explicitly terms linear in the phonon propagator D only. 
As the result one obtains the Migdal-Eliashberg theory for 

t = tc + tEP, (45) 

where 

t(l,2) = -\^"(lJ)^3G(l,2)f,(2,2;l). (46) 

( 1 , 2) = ( 1 , 2)e-i (2, 2) is the screened Coulomb interaction. The part which is due 
to the EPI has the following form 

L£p(l,2) = -y£p(l,2)fe(l,3;l)G(3,4)f,(4,2;2), (47) 

where 

y£p(l,2) = e,-i(l,I)D(l,2)e-i(2,2) 

is the screened EPI and is the electronic dielectric function. Note, L£:p(l, 2) depends 
now quadratically on the vertex function , due to the adiabatic theorem. If F^ (which 
is a functional of G) is known then the quasiparticle dynamics can be in principle deter- 
mined. In that respect the central question is: (1) how to calculate F^ - which contains 
all information on Coulomb interaction and electronic correlations? This is a difficult 
task and practically never realized in its full extent for real systems. However, this pro- 
gram is realized recently in the t-J model with the EPI interaction in the framework 
of the X-method - see below and [2]; (2) how to calculate the effective EPI potential 
Vep ~ gpp/^e^ or more precisely the coupling constant gpp and the electronic dielectric 
function gg? In absence of a better theory these quantities are usually calculated in the 
framework of the LDA band-structure theory. 



3.2. LDA calculations of X in HTSC oxides 



The LDA method considers electrons in the ground state (there is a generalization to 
finite r), whose energy can be calculated by knowing the spectrum {e;^} of the Kohn- 
Sham (Schrodinger like) equation 

[|^ + y,(r)]vAfc(r) = e^VA^(r), (48) 

which depends on the effective one-particle potential 

Y^ir) = Vei{r)+VH{r)+Vxc{r). (49) 

Here. Vei is the electron-lattice potential, Vh is the Hartree term and Vxc describes 
exchange-correlation effects - see [2]. Because the EPI depends on the excited states 
(above the ground state) of the system this means, that in principle the LDA method can 
not describe it - see [2]. However, by using an analogy with the microscopic Migdal- 
Eliashberg theory one can define the EPI coupling constant g(^'s) = gTc/s also in the 
LDA theory - see [2]. It reads 

= E«'(l'.k') = (V. I I (50) 

where n means summation over the lattice sites, a = x,y,z and the wave function 1//^/ is 
the solution of the Kohn-Sham equation. Formally one has SV^/SR,, = TiDA^e^^yei- 

Even in such a simplified approach it is difficult to calculate g'g^^^^ — g^J^ + g'an"'^ 
because it contains the short-range (local) coupling 

8a,n ~ 8a,n (k, K ) ~ (7/^ | ra \ Yvrrj) (51) 
with A/ = 1, and the long-range coupling 

^■(k,kO ~ (7/, I (R«-RO)a I W (52) 

with A/ = 0. In most calculations the local term g^^^ is calculated only, which is justi- 
fied in simple metals only but not in the HTSC oxides. In HTSC oxides the latter gives 
very small EPI coupling X^^^ ~ 0. 1, which is apparently much smaller than the exper- 
imental value A > 1 giving rise to the pessimistically small T^ [113]. The small X^^^ 
was also one of the reasons for abandoning the EPI as pairing mechanism in HTSC ox- 
ides. At the beginning of the HTSC era the electron-phonon spectral function a^F{(o) 
for the case La2-xSrxCu04 was calculated in [114] by using the first-principles band 
structure calculations and the nonorthogonal tight-binding theory of lattice dynamics. It 
was obtained A = 2.6 and for assumed jU* = 0. 13 gave Tc = 36 K. However, these calcu- 
lations predict a lattice instability for the oxygen breathing mode near Lai ssSro i5Cu04 
that is never observed. Moreover, the same method was applied to YBa2Cu30j in [1 15] 
where it was found A = 0.5 which leads at most to Tc = (19 — 30) K. In fact the calcu- 
lations in [114], [115] do not take into account the Madelung coupling (i.e. neglect the 
matrix elements with A/ = 0). 



However, because of the weak screening of the ionic (long-range) Madelung coupling 
in HTSC oxides - especially for vibrations along the c-axis, it is necessary to include 
the nonlocal term g^^'n"'^ ■ This goal was achieved in the LDA approach by the Pickett's 
group [51], where the EPI coupling for Lai^xMxCuO^ is calculated in the frozen-in 
phonon (FIP) method. They have obtained A = 1.37 and (Oiog ~ 400 K and for /i* = 0. 1 
one has Tc ^ 49 K (Tc ^ a)[ogexp{-l/[(A/(l -l-A)) - jU*]}). For more details see Ref. 
[2] and references therein. We point out, that some calculations which are based on the 
tight-binding parametrization of the band structure in YBa2Cu^0'] gave rather large EPI 
coupling X ^2 and Tc = 90 K 

Recently, a new linear-response full-potential linear-muffin-tin-orbital (LR-LMTO) 
method for the calculation of A^"* was invented in [116]. It is very efficient in ex- 
plaining the physics of elemental metals, like Al,Cu,Mo,Nb, Pb, Pd, Ta and V with dis- 
agreements by only 10 — 30% of theoretical and experimental results (obtained from 
tunnelling and resistivity measurements) for the EPI coupling constants A and Xtr- How- 
ever, the LR-LMTO method applied to the doped HTSC oxide {Cai^xSrx)i-yCu02 for 
X ~ 0.7 and y ~ 0. 1 with Tc =110^ gives surprisingly small EPI coupling As ~ 0.4 for 
s — wave pairing and A^ < 0.3 for d — wave pairing [1 17]. Although this finding, that Xd 
is of the similar magnitude as A^ {X^ A^), is interesting and encouraging it seems that 
this method misses some ingredients of the ionic structure of the layered structure [52], 
[53]. 

We point out, that the model calculations which take into account the long-range 
ionic Madelung potential appropriately [47], [48], [50] gave also rather large coupling 
constant A ~ 2, what additionally hints to the importance of the long-range forces in the 
EPI. 

Since in HTSC oxides the plasma frequency along the c-axis, (0^^ is of the order 
(or even less) of some characteristic c-axis vibration mode, it is necessary to include 
the nonadiabatic effects in the EPI coupling constant, i.e. its frequency dependence 
ga,n ~ g^/£cci(o)- This non-adiabaticity is partly accounted for in the Falter group [52], 
[53] by calculating the electronic dielectric function along the c-axis ecc(k, co) in the 
RPA approximation. The result is that ga,n is increased appreciable beyond its (well 
screened) metallic part, what gives a large increase of the EPI coupling not only in the 
phonon modes but also in the plasmon one. This question deserves much more attention 
than it was in the past. 



3.3. Lattice dynamics and EPI coupling 

The calculation of the phonon frequencies (Op/,, which are obtained from 

Do\q, COph) - fl(q, cOph) = 0, (53) 

is in principle even more complicated problem than the calculations of the electronic 
properties. It lies on the difficulty to calculate the phonon polarization operator n - see 
more in [2]. Schematically one has 

n^{VaVe-i)Hc, (54) 



where Xc is the electronic charge susceptibility. Vg-t is the bare electron-lattice interac- 
tions. Xc is schematically given by Xc = P^e^ electronic polarization operator 
P = GTcG. As we see the phonon frequencies depends crucially on the screening prop- 
erties of electrons. The screening effects in HTSC oxides are determined by the speci- 
ficity of the metallic-ionic structure and strong electronic correlations. At present there 
is a controllable theory for the electronic properties in the t-J model [17], [18], [19], 
[20], [2] only, where these two ingredients are successfully incorporated in the theory. 
However, until now there is no controllable theory for the lattice dynamics which in- 
corporates these two ingredients, in spite the fact that the X-method (see below) offers 
well defined and procedure. There were a number of interesting attempts to calculate 
renormalization of some specific phonons [54], [55], [56], such as for instance of the 
half-breathing phonon mode along the (1 , 0, 0) direction - which is strongly softened. In 
spite of some alleged theoretical confirmation of the experimental softening in YBCO 
and LASCO, none of these calculations are reliable, because none of them take into ac- 
count the screening due to strong correlations (the charge vertex and the dielectric 
function Eg) in a controllable way. That is the reason that all attempts until were unable 
to extract the reliable magnitude of the coupling constant with a specific phonon. Even 
more, by playing only with a single phonon mode, and with a particular wave- vector in 
the Brillouin zone, one can not get large EPI and large X - see [2]. The latter claim is 
confirmed by tunnelling experiments, which demonstrate that almost all phonons (for in- 
stance 39 modes in YBCO) contribute to X . No particular mode can be singled out in the 
spectral function a^F{o)) as being the only one which dominates in pairing mechanism 
in HTSC oxides. 



4. THEORY OF STRONG ELECTRONIC CORRELATIONS 

The well established fact is that strong electronic correlations are pronounced in HTSC 
oxides, at least in underdoped systems. However, the LDA theory fails to capture effects 
of strong correlations by treating they as a local perturbation. This is, as we shall see 
later, an unrealistic approximation in HTSC oxides, where strong correlations introduce 
non-locality. The shortcoming of the LDA is that in the half-filling case (with n = 1 and 
one particle per lattice site) it predicts metallic state missing the existence of the Mott 
insulating state. In the latter, particles are localized at lattice sites independent of the 
(non)existence of the AF order and the localization is due to the large Coulomb repul- 
sion ?7 at a given lattice site, i.e. U where W is the band width. Some properties 
in the metallic state can not be described by the simple canonical Landau-Fermi liquid 
concept. For instance, recent ARPES photoemission measurements [122] on the hole 
doped samples show a well defined Fermi surface in the one-particle energy spectrum, 
which contains 1 — 5 electrons in the Fermi volume (5 is the hole concentration), but 
the band width is (2 — 3) times smaller than the LDA band structure calculations predict. 
The latter is consistent with the Luttinger theorem as well as with the LDA band struc- 
ture calculations. However, experimental data on the dynamical conductivity (spectral 
weight of the Drude peak). Hall measurements etc. indicate that in transport proper- 
ties a low density of hole-like charge carriers (which is proportional to 5) participates 



predominantly. These carriers suffer strong scattering and their inverse lifetime is pro- 
portional to the temperature (at T > Tc) as we discussed earlier. It is worth of mentioning 
here that the local moments on the Cu sites, which are localized in the parent AF com- 
pound, are counted as part of the Fermi surface area when the system is doped by small 
concentration of holes in the metallic state. The latter fact gives rise to a large Fermi 
surface which scales with the number (per site) of electrons 1 — 5. At the same time 
the conductivity sum-rule is proportional to the number of doped holes 5, instead of 
1 — 5 as in the canonical Landau-Fermi liquid. These two properties tell us that we deal 
with a correlated state, and the latter must be due to the specific electronic structure of 
HTSC oxides (cuprates). The common ingredient of all cuprates is the presence of the 
Cu atoms. In order to account for the absence of Cm^"*" ionic configuration (the charge 
transfer Cu^^ Cu^^ costs large energy f/ ~ 10 eV , i.e. the occupation of the Cu site 
with two holes with opposite spins is unfavorable) P. W. Anderson [61] proposed the 
Hubbard model as the basic model for quasiparticle properties in these compounds. For 
some parameter values it can be derived from the (minimal) microscopic three-band 
model. Besides the hopping tpd between the J-orbitals of Cu and p-orbitals of O ions 
(as well as tpp) - the Emery model [123], it includes also the strong Coulomb interaction 
Uqu on the Cu ions as well as interaction between p- and d-electrons. The main two pa- 
rameters are Ucu ~ (6 — 10) eV and the charge transfer energy Ap^ = — ej^ ~ (2.5 — 4) 
eV, where e^, are energies of the d- and p-level, respectively. In HTSC oxides the case 
Ucu » ^pd is realized, i.e. they belong to the class of charge transfer materials. This 
allows us to project the complicated three-band Hamiltonian onto the low-energy sector, 
and to obtain an effective single-band Hubbard Hamiltonian with an effective hopping 
parameter t and the effective repulsion U Ap^. It turns out that the case U >> ? is 
realized, since Apj ^ t = t^^/Apd. The effective and minimal Hamiltonian which de- 
scribes the low-energy physics of HTSC oxides comprises also the long-range Coulomb 
interaction Vc and the EPI Vepi - see [7], [2] 

^=-Jl tiflcjCjc + L m^uii + Vc + Vepi. (55) 

i,j,a i 

The effective repulsion U ^4 eV has its origin in the charge-transfer gap of the three- 
band model, while the nearest neighbor and next-nearest neighbor hopping t and t' , 
respectively are estimated to be ? = 0.3 — 0.5 eV and t'/ 1 equal —0. 15 in La compounds 
and —0.45 YBCO. Since (i7/?) 3> 1 the above Hamiltonian is again in the regime of 
strong electronic correlations, where the doubly occupancy of a given lattice site is 
strongly suppressed, i.e. -C 1. The latter restricts charge fluctuations of electrons 

(holes) on a given lattice site are allowed, since = 0, 1 is allowed only, while processes 
with ni — 2 are (practically) forbidden. Note, that in (standard) weakly correlated metals 
all charge fluctuation processes (n, = 0, 1 , 2) are allowed, since U <^W in these systems. 
From the Hamiltonian in Eq.(55), which is the 2D model for the low-energy physics in 
the Cu02 plane, comes out that in the undoped system there is one particle per lattice 
- the so called half-filled case (in the band language) with (n,) = 1. It is an insulator 
because of large U and even antiferromagnetic insulator at T = ^. The effective 
exchange interaction (with J = At'^ /U) between spins is Heisenberg-like. By doping 
the system by holes (with the hole concentration 5{< 1)), means that particles are taken 



out from the system in which case there is on the average (n,) = 1 — 5 particles per 
lattice site. Above some (small) critical doping 5c ~ 0.01 the AF order is destroyed and 
the system is strongly correlated metal. For some optimal doping 5op(~ 0. 1) the system 
is metallic with the large Fermi surface and can exhibit even high-Tc superconductivity 
in the presence of the EPI, as it will demonstrated below. The latter interaction and its 
interplay with strong correlations is the central subject in the following sections. 



4.1. X-method for strongly correlated systems 

Since U »t one can put with good accuracy 11^°°, i.e. the system is in the strongly 
correlated regime where the doubly occupancy = 2 is excluded. One of the ways 
to cope with such strong correlations is to introduce the (fermionic like) creation and 
annihilation operators (Zf ° and XP^ = (Xf 0)t) 

Z/^o = 4(l-n,,_^), (56) 

which respect the condition o- + < 1 on each lattice site. The latter means that 
there is no more than one electron (hole) at a lattice site, i.e. the doubly occupancy is 
forbidden. The bosonic like operators 

(with o\ ^ 02) create a spin fluctuation at the / — th site. Here, the spin projection 
parameter cr =t, J, and — a =i, t and the operator X^^ has the meaning of the electron 
(hole) number on the i-th site. In the following we shall use the convention that when 
^i^^ I 1) — 1 I 1) there is a fermionic particle ("electron") on the /-th site, while for 
^f*^ I 0) = I 0) sit^ is empty, i.e. there is a hole on it. It is useful to introduce the 
hole number operator 

X«o = xO^ZfO (58) 

at a given lattice site, i.e. if Xf^ | 0) = 1 | 0) the z-th site is empty - there is one hole on 
it, while for Xf^ \ 1 ) = | 1 ) it is occupied by an "electron" and there is no hole. They 
fulfill the non-canonical commutation relations 



xf,x]\ = 5,j[5,pX-'±5,,xf 



(59) 



Here, ct, j8 , y, A = 0, a and 5,y is the Kronecker symbol. The (anti)commutation relations 
in Eq.(59) are rather different from the canonical Fermi and Bose (anti)commutation 
relations. 

Since U ^ °° the doubly occupancy is excluded, i.e. | V)(= Xf^ iTi)) =0, and 
by construction the X operators satisfy the local constraint (the completeness relation) 

Cx(0^#°+i;#'^=l. (60) 

(T=l 



This condition tells us that at a given lattice site there is either one hole (Xj^^ \ hole) = 
1 I hole)) ore one electron (X-^*^ | elec) = 1 | elec)). Note, if Eq.(60) is obeyed then both 
commutation and anticommutation relations hold also in Eq.(59) at the same lattice site, 
which is due to the projection properties of the Hubbard operators X^^X'^^ = d^yX"^. 

For further purposes, i.e. for the study of low-energy excitations in a controllable way, 
the number of spin projections is generalized to be , i.e. o = 1,2,...N. This way of 
generalization was very useful in describing heavy fermion physics, where for some 
Ce compounds means the number of projections of the total angular momentum, for 
instance when j = 5/2 then A^ = 2j -|-1 = 6 (N 1). For some Yb compounds one 
has j = 7/2, i.e. N = 2j + \ = S (N ^ I). By projecting out the doubly occupied (high 
energy) states from the Hamiltonian in Eq.(55) one obtains the generalized t — J model. 
The details of this derivation are given in Appendix and here we give the final expression 
for the Hamiltonian which excludes the doubly occupancy 

H = H,+Hj = -l^ tijXr'xf^ +l^j.j^SfSj-^-hihj)+H3. (61) 

The first term describes the hopping of the "electron" by taking into account that the 
doubly occupancy of sites are excluded. The second term describes the Heisenberg- 
like exchange energy of almost localized "electrons". H3 contains three-sites hopping 
which is usually omitted believing it is not important. For effects related to charge 
fluctuation processes it is plausible to omit it, while for spin-fluctuation processes it 
may be questionable approximation. The spin and number operators can be expressed 
via the Hubbard operators [125] 

S = xf'°(d)^,^,X°*^;n,- = zr (62) 

where summation over bar indices is understood. 

The basic idea behind the X-method is that the Dyson's equation for the electron 
Green's function can be effectively obtained by introducing external potentials (sources) 
u<^i<^2(^iy fiiQ source Hamiltonian Hg is used in the form 

[H,dZ= f £ M^"^2(1)X^1^2(1)C/1 = M*1*2(1)X*1*2(1), (63) 
J J <7l,<72 

where 1 = (1, t) and J{..)dl = /(-O^^Li ^nd T is the Matsubara time. Here and 
in the following, integration over bar variables(l,2..) and a summation over bar spin 
variables(d'..) is understood. The sources u^^^^{l) are useful in generating higher corre- 
lation functions entering the self-energy. The electronic Green's function is defined by 
[2] (f is the time-ordering operator) 

, -(f (5X0^1(1)1^20(2))) 
{1 i) 

where S = fexp{— J Hs{l)dl} and the corresponding Dyson's equation reads 

Gq^^' ( 1 , 2) - *2 ( 1 , 2) j G^2ct2 (^2, 2) = ''H 1 ) ^ ( 1 - 2) , (65) 



G-;'^>^2(l,2) = (-^)5^i^25(l-2)-M^i*^2(i)5(l-2) (66) 
The so called Hubbard spectral weight is given by 

(2«^i^2(i) = 5''i^2^x00(i)) + (1^1^2(1)). (67) 

Eg' '^^(1,2) is a functional of the Green's function G'^^'^^{1,2). The latter describes the 
composite (correlated) particle (in the language of the 55 theory it describes the com- 
bined "spinon + holon"). For further analysis it is useful to introduce the quasiparti- 
cle Green's g^i*^2 (something analogous to the "spinon" Green's function in the SB ap- 
proach) and the vertex functions Jalal ( 1 7 2; 3) , respectively 

^<7i<72(i^2) = G*''^(1,2)!2"^'^2*'H2) (68) 

X -l,C7i(T2('1 2) 

ygg(l,2;3) = - "^,..,,,^,r\ (69) 



respectively. Note, that 703^4 (1,2; 3) are the three-point vertex function, which also 
renormalizes the ionic EPI coupling constant - as we shall see below. g^^^^{l,2) is the 
solution of the equation 

^o^i'""'*" ( 1 , 2) - Eji ^2 ( 1 , 2) 1 ^^2^2 (2, 2) = ^2 5 ( 1 _ 2) , (70) 

where Eg ' '^^{1,2) depends on the "quasiparticle" Green's function g'^^ ^^(1,2). Note, that 
in Eq.(70) the Hubbard spectral weight Q^*^ disappears from the right hand side. 

Since in the following we study nonmagnetic (paramagnetic) normal state one has 
EJ^(1,2) = Eg(l,2) for C7 = 1, ..A^ , as well as Q''''{ = Q). In [17], [18], [19], [20] it is 
shown that in that case Eg can be expressed via two vertex functions - the charge vertex 



and the spin vertex 
i.e. Eg is given by 



7,(1,2;3) = 0^?(1,2;3) (71) 
7,(l,2;3) = o|^(l,2;3), (72) 



Eg(l -2) = J^^Qil) + 5il-2)^^<X-<^(l) > 
'°^^"^^g(T-2)7c(2,2;l)+ 



N 

^^^^g{l -2)7.(2,2;3) +E2(1 -2), (73) 

where ^2(1, 2, 3) = 5(1 - 2)?o(l - 3) - 5(1 - 3)/o(l - 2). The notation to{l - 2) (and 
7o(l — 2)) means to(l — 2) = toj^j25{Ti — T2). The first two terms in Eq.(73) represent 
an effective kinetic energy of quasiparticles in the lower Hubbard band. As we shall 



see below they give rise to the band narrowing and to the shift of the band center, re- 
spectively. The third and fourth terms describe the kinematic and dynamic interaction of 
quasiparticles with charge and spin fluctuations, respectively, while the very important 
term proportional to Eg (1,2) takes into account the counterflow of surrounding quasi- 
particles which takes place in order to respect the local constraint (absence of doubly 
occupancy). It reads 



V n o^_ ^o(l-l)g(l-2) . dQ{2) 5Q--i2) , 



L depends on the vertex functions 7c (1,2; 3) and 7^(1, 2; 3) and it does not contain a 
small expansion parameter, like the interaction energy in weakly interacting systems, 
because the hopping parameter t describes at the same time the kinetic energy and kine- 
matic interaction of quasiparticles. This means that there is no controllable perturbation 
technique, due to the lack of small parameter. There are various decoupling procedures 
and mean-field like techniques - the path integral method, or 1 /N expansion in various 
slave-boson approaches [2]. 

What is the advantage of the X-method expressed by Eq.(68-73). It turns out that 
it allows to formulate a controllable 1/A^ expansion for Eg by including also the EPI 
[17], [18], [19], [20] - see below. For that purpose it is necessary to generalize the local 
constraint condition 

Cx,(/)^#0+i;xr = y, (74) 

(T=l ^ 

where N/2 replace the unity in Eq.(60). It is apparent from Eq.(60) that forN = 2 it 
coincides with Eq.(60) and has the meaning that maximally half of all spin states at a 
given lattice site can be occupied. 

The spectral function A(k, co) — — ImG(k, co)/n must obey the generalized Hubbard 
sum rule which respects the new local constraint in Eq.(74). 



/ 



d(0A{k, CO) = ^^^^ (75) 



The N > 2 generalization of the local constraint allows us to make a controllable 

l/N expansion of the self-energy with respect to the small quantity 1/A^ (when A'^ ^ 1). 
Physically this procedure means that we select a class of diagrams in the self-energy 
and response functions which might be important in some parameter regime. By careful 
inspection of Eq.(73) one concludes that for large there is 1 /N expansion for various 
quantities - for instance 

g^gO + ^ + -;Q^Nqo + Qi + ..., (76) 

Eg = Eo + ^ + ...;7c = 7cO + |^ + -, (77) 
7,(l,2;3)=iV5(l -2)5(1 -3) + 7,i + ... 



As a result of this expansion one obtains Eq and go in leading 0{\)-order - see details 
in [2], [17], [18], [19], [20] (Eg in Eq.(73) is of 0(1/N) order) 



^o(k, CO) 



Gojk, a) 
Go 



(0-[eo(k)-M]' 



1 



(78) 



where the quasiparticle energy eo(k) = — ^o?o(k) — Lp-^o(k + p)nF(p), and the level 
shift £c = Lp^o(p)"f(p)- Here, ?o(k) and /o(k) are Fourier transforms of tojj and Jo ij, 
respectively. For t' = one has ?o(k) = 2to{cos kj^ + cos ky) and io(k) = 2jQ{coskx + 



Let us summarize the main results which were obtained by the X-method in leading 
C>(l)-order and compare these results with corresponding results of the 5fi-method 
[17], [18], [19], [20]: (1) In the 0(1) order the Green's function go(k,w) describes 
the coherent motion of quasiparticles whose contribution to the total spectral weight of 
the Green's function Go(k, (o) is Qq = N5/2. Note, Go(k, co) = 2o^o(k, (o) in leading 
order. The dispersion of the quasiparticle energy is dominated by the exchange parameter 
if 7o > ^^0- In the case when Jq — there is a band narrowing by lowering the hole 
doping 5, where the band width is proportional to the hole concentration 5, i.e. W = 
z - 5 ■ tQ. (2) The X-method respects the local constraint at each lattice site and in 
each step of calculations. (3) In the important paper [126] - which is based on the 
theory elaborated in [17], [18], [19], [20], it is shown that in the superconducting 
state the anomalous self-energy (which is of 0(l/A^)-order in the 1/A^ expansion) of 
the X- and 55-methods differ substantially. As a consequence, the 55-method [127] 
predicts false superconductivity in the t — J model (for / = 0) with large T^ (due to the 
kinematical interaction), while the X-method gives extremely small rc(pa 0) [126]. The 
reason for this discrepancy between the two methods is that calculations done in the 
S5-method miss a class of compensating diagrams, which are on the other hand taken 
automatically in the X-method. So, although the two approaches yield some similar 
results in leading 0(l)-order their implementation in the next leading 0(1/A^) -order 
make that they are different. Note, that the 1 /N expansion in the X-method is well- 
defined and transparent. (4) By explicit calculation and comparison of the two methods 
in [18] it is shown that the renormalization of the EPI coupling constant is different in the 
two approaches even in the large A''- limit - see below; (5) Very interesting behavior as a 
function of doping concentration exhibits the optical conductivity (7(0), q = 0) = <7((o) 
which scales with the doping 5. Note, the volume below the Fermi surface in the 
case of strong correlations scales with n = 1 — 5, like in the usual Fermi liquid. The 
above analysis clearly demonstrate shows difference in response functions of strongly 
correlated systems and the canonical Landau-Fermi liquid. 




(79) 



The chemical potential is obtained from the condition 



(80) 



p 



4.2. Forward scattering peak in tlie cliarge vertex 7c 



The three-point charge vertex 7c( 1,2; 3) plays important role in the renormalization 
all charge processes, such as the EPl, Coulomb scattering and the scattering on non- 
magnetic impurities. It was shown in [17], [18], [19], [20], [2] that 7^(1, 2;3) can be 
calculated exactly in the leading 0(1) order (and in all other orders) of the ? — 7 model. 
The integral equation in the C>( 1) order reads 

7(1, 2;3) = 5(1-2)5(1-3) +?(1 -2)go(l, l)go(2, 1+)7(1,2;3) 

+5(1 -2)f(l - l)go(l,2))go(3, 1)7(2,3;3 

-7(l-2)go(l,l))^o(2,2)7(l,2;3. (81) 
The analytical solution of Eq.(81) is given by [2], [20] 

7c(k,^) = l- £ ZFa{^)[^+x{q)]alXp2{q). (82) 
a=lj8=l 

where 

Xafi{q)^Y.'^a{p,q)F^{v), (83) 
P 

Gaip,q) = [^,tip + q), cos px, sin px, cos py, sin py]n{p,q), 

Fa(k) = [t{k),\,2Jocoskx,2Josinkx,2Jocosky,2Josinky], 

where U{k,q) = -g{k)g{k + q) and q = {q,iqn), qn = 2nnT, p = {p,ipm), Pm = 
TtT{2m-\- 1). Note, the frequency sum over pm in Xapiq) in Eq.(82) involves only IT 
and can easily be carried out Y.p^ ^ip,q) = ["^('^q+p) - «f('^p)]/['^p - ^q+p - iqn]- 

We stress that Yc{kjq) describes a specific screening of the charge potential due to of 
strong correlations. In the presence of perturbation (external source u) there is change 
of the band width, as well as of the local chemical potential, which comes from the sup- 
pression of doubly occupancy. The central result of the X-method is that for momenta 
k laying at (and near) the Fermi surface 7.o(k, q, co = 0) has very pronounced /orwarJ 
scattering peak at q = 0) at low doping concentration 5(<^ 1), while the backward scat- 
tering is substantially suppressed - see Fig. 17. The latter means that charge fluctuations 
are strongly suppressed (correlated) at small distances. Such a behavior of the vertex 
function means that a quasiparticle moving in the strongly correlated medium digs up 
a giant correlation hole with the radius t,ch'^a/8, where a is the lattice constant - see 
Fig. 16. 

However, in the highly doped systems with 5 > 0.1 - which corresponds to the 
overdoped HTSC oxides, the effects of strong correlations is progressively suppressed 
and the screening mechanism due to strong correlations is less effective. We stress that 
when J <t then the last term in Eq.(81) for 7co(kir,q) is unimportant. On the other 
hand both terms, the second (due to band narrowing) and the third (due band shifting) 
one, are in conjunction responsible for the development of the forward scattering peak 
at lower doping. If we omit in Eq.(81) the band shifting term (the third one) we get very 




FIGURE 16. Schematic picture of electron correlation hole and the E — P interaction for uncorrelated 
(weakly correlated) (a) and strongly correlated (b) electron. In the case (a) the electron does not perturb 
the electronic density n(r) and it interacts with the vibrations of a single atom (shaded). In the case (b) an 
electron is accompanied by a large correlation hole of size ^ ^ 1/5 (5 is doping) and it will interact with 
atoms within this zone. From [120]. 



weak foi^ard scattering peak, while omitting the band narrowing term (the first one) 
Ycoi^F, q) is practically constant in a broad region of q. 

Finally, since the real physics is characterized by N — 2 one can put the question 
- what is the reliability of the results for the quasiparticle properties obtained by the 
l/N expansion (and N ^ oo) within the X-method? First, the exact diagonalization of 
the charge correlation function A^(k, (o) in the t — J model [124] shows clearly that 
the low-energy charge scattering processes at large momenta | k |Ri 2kp are strongly 
suppressed compared to the small transferred momenta (| k |<^ Ikp). These calculations 
confirm unambiguously the results obtained by the X-method in [17], [18], [19] on 
the suppression of the backward scattering in the vertex. Second, very recent Monte 
Carlo (numerical) calculations in the Hubbard model with finite U [21] show clear 
development of the forward scattering peak in Yci^F, q) by increasing U, thus confirming 
the theoretical predictions in [17], [18], [19]. 
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FIGURE 17. Zero-frequency vertex function 7(kf , q) of the f — 7 model as a function of the momentum 
aq with q = (q, q) for three different doping 5 - from [19]. 



5. RENORMALIZATION OF THE EPI BY STRONG 

CORRELATIONS 

In preceding Sections arguments we argued that because Xtr <C A the standard Migdal- 
Eliashberg EPI theory must be corrected in order to take into account screening prop- 
erties of strongly correlated system. This analysis is done in the framework of the X- 
method in a series of papers [17], [18], [19] which we briefly discuss below. The renor- 
malization of the EPI by strong correlations has been studied also by the S'fi-method 
[128], [121], [129], [130] by the 1/A^ expansion in the partition function, or by using 
the mean-field approach [131] where also the non-Migdal correction due to the EPI is 
considered. We stress that at present there are no systematic and controllable calcula- 
tions within the S5-method for the EPI. From that point of view the Z-method is of 
indispensable value. 



5.1. The forward scattering peak in the EPI 



The minimal model Hamiltonian for the HTSC oxides contains besides thet — J terms 
also the EPI, i.e. ^ = Htj + Hep where Hep = H'£^ + H^p 

^ = - I t^jXr'X^" + LM^i ■ Sy - \ntnj )+ 

+ L ^ajXr + Hph + Hep + Vlc, (84) 

i,(J 

where the ionic contribution to the EPI is 

Hk7 = i4>Kxr - (xn) • (85) 

(,<J 

Here, <!>, (given by Eq.(40)) describes the change of the atomic energy e° • due to the 
long-range Madelung energy, where L and Jc enumerate unit lattice vectors and atoms 
in the unit cell, respectively. Z^^ is the effective charge of an ion at the site Lk. Note, 
in Eq.(85) we do not assume small displacement u, and the following analysis holds 
in principle also for an anharmonic EPI. The term proportional to {X^'^) in Eq.(85) is 
introduced in order to have (4>;) = in the equilibrium state. Note, that there is also 
covalent contribution to the EPI in Eq.(84) due to the change of the hopping (?) and 
exchange energy (7) by the ion displacements 



+E 



(u,-Uj)SrSj. (86) 



The treatment of the first term is similar to the Madelung term in Eq.(85) although 
the equation for \hQ four-point vertex function 7c(l,2;3,4) is different than Eq.(86). 
We stress, that the X-method has advantage also in the treatment of the covalent term, 
because it peaks up straightforwardly all important contributions in jc, due to strong 
correlations. On the other side the corresponding treatment by the SB method is com- 
plicated and not well defined, giving sometimes wrong results. For instance, in [121] 
several terms in the vertex equation are omitted leading to incorrect results for the co- 
valent EPI coupling. We stress, that the covalent part contributes approximately 20 — 30 
K to the critical temperature in HTSC oxides as the band structure calculations in [1 14], 
[115] have shown (partly discussed in Section 3.). Its renormalization by strong corre- 
lations will be studied elsewhere. The second term in H^p is due to the change of the 
exchange energy by phonon vibrations. Since it is second order with respect to tij it is 
much smaller than the first covalent term and accordingly contributes very little to the 
total EPI. 



After technically lengthy calculations, which are performed in [17], [18] the expres- 
sion for the ionic part of the EPI (frequency-dependent) part of the self-energy reads 

L|5:''^(l,2) = -y£p(l -2)7,(1, 3;l)go(3-4)r,(4,2;2), (87) 

where analogously to Eq.(47) one has VEpil - 2) = £-\l - l)V^p(l - 2)e-^{2 - 2). 
The propagator of thebareEPI Vgp(l — 2) = — (r4>(l)<i>(2)) comprises in principle also 
the anharmonic contribution. From Eq.(87) it is seen that in strongly correlated systems 
the ionic part of the EPI is proportional to the square of the three-point charge vertex 
7c(l , 2; 3) (due to correlations). The self-energy is given by 

poo 

r^p {Ko))^ da{a^F{k,k',a))i,,R{co,a), (88) 
Jo 

where R{a), O.) is given in [2], [17], [18]. The (momentum-dependent) Eliashberg spec- 
tral function is defined by 

a^F{ky,co) =A^.c(0)£ I geff{Kk-k',v) |2 X 

V 

x5((0-(0v(k-k'))T^(k,k-k'). (89) 

is the Bose distribution function and Xj/ is di-gamma function, while ge//(k,p, v) 
is the EPI coupling constant for the V-the mode, where the renormalization by long- 
range Coulomb interaction is included, i.e. ge//(k,p, v) = ^(k,p, v)/ee(p). Nsc{0) is 
the density of states renormalized by strong correlations where Nsc{0) = NQ{0)/qo and 
qo — 5/2mthet — t' model (7 = 0). In thet — J model NsdO) has another form which 
does not diverge for 5 — > but one has Nsc{0) (~ l/7o) > No{0), where the bare density 
of states A'o(O) is calculated, for instance by the LDA scheme. 



5.2. Pairing and transport EPI coupling constants 

Depending on the symmetry of the superconducting order parameter A(k, (o) (s—, 
d — wave pairing) various averages (over the Fermi surface) of a^F{k,k\ (o) enter the 
Eliashberg equations. Assuming that the superconducting order parameter transforms 
according to the representation (/ = 1 , 3, 5) of the point group €4^ of the square lattice 
(in the Cu02 planes) the appropriate symmetry-projected spectral function is given by 

a^Fi{k,^,CD) = ^ J: I geff{k,k-Tjk',v) |2 X 

xd{co-(Oy{k-Tjk')) I rc{k,k-Tjk') l^DiU). (90) 

k and k' are momenta on the Fermi line in the irreducible Brillouin zone which is 1/8 

of the total Brillouin zone. Tj , j = 1,..8, denotes the eight point-group transforma- 
tions forming the symmetry group of the square lattice. This group has five irreducible 



representations which we distinguish by the label i = 1, 2, ...5. In the following the rep- 
resentations i — 1 and i = 3 , which correspond to the s— and d — wave symmetry of 
the full rotation group, respectively, will be of importance. Di{j) is the representation 
matrix of the j — th transformation for the representation /. By assuming that the su- 
perconducting order parameter A(k, to) does not vary much in the irreducible Brillouin 
zone one can average over k and k' in the Brillouin zone. For each symmetry one obtains 
the corresponding spectral function a^Fi{o}) 

a'Fi{(o) = {{a'Fi{kMM\h (91) 

which (in the first approximation determines) the transition temperature for the order 
parameter with the symmetry F;. In the case i = 3 the electron-phonon spectral function 
a^F2,{(o) in the d-channel is responsible for d — wave superconductivity represented by 
the irreducible representation (or sometimes labelled as B\g). 

Performing similar calculations (as above) for the phonon-limited resistivity one finds 
that the latter is related to the transport spectral function a?'Ftr{(o) which is given by 

((«^F(k,k^a))[v(k)-v(kO]2)^)^, 

2((v2(k))k)k' 

v(k) is the Fermi velocity. The effect of strong correlations on the EPI was discussed 
in [17] and more extensively in [18] within the model where the phonon frequencies 

(0(k — k ) and ^e//(k,p,A) are weakly momentum dependent - due to the long-range 
screening (RPA). In order to illustrate the effect of strong correlations on a^Fi{o}) we 
consider the latter functions at zero frequency {(O = 0) which are then reduced to the (so 
called) "enhancement" functions 

= \^ E ((I rc(k,k-7}-k') WhDiU) (93) 

Note, in the case 7 = one has A^,c(0)/A^o(0) = ^ where qo is related to the doping 
concentration, i.e. qQ = 8 /I. Similarly, the correlation effects in the resistivity p{T){r^ 
Atr) renormalize the transport coupling constant Atr) 

. _ NscjO) ((I YciU-Tjk') p [v(k) -y{k')]\)y,, 

No{0) 2((v2(k))k)k. ^ ^ 

Note, that for quasiparticles with the isotropic band the absence of correlations implies 
that Ai = A;^ = 1 , A, = for / > 1 . 

The averages in Ai,A3 and A^^ were performed numerically in [18] by using the 
realistic anisotropic band dispersion in thet — t' — J model and the corresponding charge 
vertex. The results for Ai,A3 and A,^ are shown in Fig. 18 as functions of doping 
concentration in the t and t — t' and t — t' — J models, respectively. The three curves 
are multiplied with a common factor so that A] approaches 1 in the empty-band limit 
5^1, when strong correlations are absent. Note, that Tc in the weak coupling limit and 
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FIGURE 18. (a) - Enhancements Ai and A3 and 5 • Afr as a function of doping 5 for t' =0 and 7 = 
- from [18]. (b) Enhancements Ai and A3 and 5 • Atr as a function of doping 5 for t' = —0.05 and 7 = 0- 
from [18]. 



in the / — th channel scales like 



where /if is the Coulomb pseudopotential in the i-th channel and {(o) averaged phonon 
frequency. 

Several interesting results, which are seen in Fig. 18, should be stressed. 

First, in the empty-band limit 5 ^ 1 the J — wave coupling constant A3 is much 
smaller than the s — wave coupling constant Ai, i.e. A3 ^ Ai. Furthermore, the totally 
symmetric function Ai decreases with decreasing doping. 

Second, in both models Ai and A3 meet each other (note Ai > A3 for all 5) at some 
small doping 5 ^0.1 —0.2 where Ai A3 but still Ai > A3. By taking into account a 
residual Coulomb repulsion of quasiparticles with 11*1 <^ ji* one gets that the s — wave 
superconductivity (which is governed by the coupling constant Ai) is suppressed, while 
the d — wave superconductivity (governed by A3) is only weakly affected. In that case 
the d — wave superconductivity due to the EPI becomes more stable than the s — wave 

superconductivity at sufficiently small doping 5, i.e. > Tc''\ Experimentally, this 
occurs in underdoped, optimally doped an overdoped HTSC oxides [34]. This transition 
between s- and d — wave superconductivity is triggered by electronic correlations 
because in the calculations it is assumed that the bare EPI coupling is momentum 
independent, i.e. the bare coupling constant contains the s — wave symmetry only. 

Third, the calculations are performed in the adiabatic approximation, where Ai is less 
and Afr much more suppressed by strong correlations. In the nonadiabatic regime (o > 
P • Vf(p) i-e. for (0< p • Vp-(p) the vertex function grows by decreasing q finally reaching 
7t (kF,p = 0,0)) = 1. Due to the latter effect the enhancement function 7^(k/7,p,(o)/^o 
may be substantially larger compared to the adiabatic. This means that different phonons 




((0)exp(-l/(A,-Ai;), 



(95) 



will be differently affected by strong correlations. For a given frequency the coupling to 
phonons with momenta p < pc — (o/vp will be enhanced, while the coupling to those 
with p > pc = (o/vp '^^ substantially reduced due to the suppression of the backward 
scattering by strong correlations. 

Fourth, the transport coupling constant Atr (not properly normalized in Fig. 18 - 
see correction in [19]) is reduced in the presence of strong correlations, especially for 
lower doping where Atr < A/3. This is very important result because it resolves the 
experimental puzzle that Xtr (which enters resistivity p (T) ~ XtrT) is much smaller than 
the coupling constant A (which enters the self-energy E and Tc), i.e. why Af^ << A. The 
answer lies in strong correlations which causes the forward scattering peak in charge 
scattering processes - the FSP theory. 

As we already said, Monte Carlo (numerical) calculations if the Hubbard model at 
finite U - performed by Scalapino Group [21], show that the forward scattering peak in 
the EPI coupling constant (and the charge vertex) develops by increasing U. The latter 
effect is more pronounced at lower doping. The similar (to Monte Carlo) results were 
obtained quite recently in [22] in the framework of the Riickenstein-Kotliar (four slave- 
boson) model. These numerical results prove the correctness of the EPI theory based on 
the X-method. 

We stress that contrary to the X-method, where the systematic 1 /N calculations of 
the EPI self -energy is uniquely done, this is still a problem for the SB (Barnes) method 
where the 1 /N expansion of the partition function Z{T,ix) is usually performed [121]. 
The existing expression (in the literature) for the vertex function in the SB method is 
different than that in the X-method [18], [18]. It seems that such a not well-controlled 
procedure omits a class of diagrams giving inadequate behavior of the coupling constant 
A as a function of doping. Additionally, the vertex function in the SB approach is peaked 
not at ^ = but at some finite ^max> where ^max only for doping 5 - see [23]. 



6. FSP THEORY AND NOVEL EFFECTS 

There are a number of effects which are predicted by the FSP theory. We have already 
explained the effects of the forward scattering peak on the EPI. In Section 2. the shift- 
puzzle in ARPES was also explained by the FSP theory (model). We discuss briefly 
some other predictions of the FSP theory containing parts not comprised in [2]. 



6.1. Nonmagnetic impurities and robustness of d-wave pairing 

In the presence of strong correlations the impurity potential is also renormalized and 
the effective potential in the Born approximation is given by M^(q) = 7c (^f, q)Mo(l)' 
where uo{q) is the single impurity scattering potential in the absence of strong corre- 
lations [142]. Since the charge vertex 7t(pF,q) is peaked at q = the potential u{q) 
is also peaked at q = 0. This means that the scattering amplitude contains not only the 
s-channel (as usually assumed in studying impurity effects in HTSC oxides), but also 
the d-channel, etc. Based on this property the FSP theory succeeded in explaining some 



experimental facts, such as: (i) the suppression of the residual resistivity p, [17], [18]. It 
is observed in the optimally doped YBCO, where the resistivity p(r) at T = ^ has a 
rather small value < 10 ilD.cm.; (ii) the robustness of <i — wave pairing [142]. The previ- 
ous theories [140], which assume M(q) = const, i.e. the s-wave scattering chaimel only, 

predict that Tc(Pi^c) = at much smaller residual resistivity p^"] ~ 50 {iQ-cm, while the 

experimental range is 200 iiQ.cm < p^^^ < 1500 ilQ.cm [141]. The latter experimental 
fact means that d-wave pairing in HTSC is much more robust than the standard theory 
predicts, and it is one of the smoking gun experiments in testing the concept of the for- 
ward scattering peak in the charge scattering potential. It is worth of mentioning that in 
a number of papers the pair-breaking effect of non-magnetic impurities in HTSC was 
analyzed in terms of the impurity concentration n,, i.e. the dependence Tc{ni). However, 
Hi is not the parameter which governs this pair-breaking effect. The more appropriate 
parameter for discussing the robustness of d-wave pairing is the impurity scattering am- 
plitude r(0, d'), which can be related to the measured residual resistivity pi which leads 
to the dependence of Tc{pi). The robustness of d-wave pairing in HTSC can be revealed 
only by studying the experimental curve Tc{pi), what has been first recognized experi- 
mentally in [141] and theoretically in [17], [142]. 

The theory of the robustness of d-wave pairing in HTSC was elaborated first in 
[142], where the FSP theory [17], [18] is applied to this problem. We shall not go 
into details - which are given in [142], [2], but we give here a general formula for 
the Tc{pi) dependence in anisotropic (including unconventional) superconductors, only. 
We assume that in anisotropic superconductivity the superconducting order parameter 
has the form A(e) = AoF(e) and generally one has {Y{0)) ^ ({Y*{e))Y{e)) = 1) 
where the momentum dependent impurity scattering amplitude is r(0, 6') = (0, 6') + 

raYd{e)Yd{e') + 

ln|^ = ^(^) - ^(^ + (1 - J8)x) - (F(0))2[^(l) i)]. (96) 

Here, X = Ts/4nTc, /3 = r^/Fs and {Y{6)) means an averaging over the Fermi surface. 
Note, Eq.(96) holds independently of the scattering strength F^,Frf, i.e. it holds in the 
Bom as well as in the unitary limit. The residual resistivity p/ can be related to the 
transport scattering rate by p, = 47trtr/ (O^i, while the s-wave amplitude is related to Ttr 
by Ts = pTtr- The parameter p > \ can be obtained from the microscopic model (for 
instance in the t-J model p fa 2 — 3) or can be treated as a fitting parameter - see more in 
[2]. In the case of an unconventional pairing one has (y (0)) = and the last term drops. 
For the s-scattering only (Fj = 0) one has j8 = and Tc{pi) should be suppressed very 
strongly contrary to the experimental results [141] - see Fig. 20. 

The FSP theory of the impurity scattering in the t-J model [142] gives that the s- 
chaimel and d-chaimel almost equally contribute to the impurity scattering amplitude, 
since j8 ^ 0.75 — 0.85 for doping 5 0.2. The dependence of j8 (5) is calculated 

for the t-J model - see Fig. 19. 

Since the d-channel in scattering is not detrimental for d-wave pairing the FSP theory 

(FSP) (s) 

predicts that Tc{pi) vanishes at much larger p,^c> i-e. p-^ » p-J, what is in the good 
agreement with experiments - as shown in Fig. 20. 
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FIGURE 19. The anisotropy scattering parameter /3 as a function of doping 5 in the t-J model. From 
[142]. 

6.2. Transport properties and superconductivity 

The EPI was studied in the past in the extreme limit of the forward scattering peak 
in the Einstein model with the phonon frequency [119], where in leading order 
the spectral function is singular, i.e. a^F(k,k',(o) ~ 5(k — k')5((0 — H). Numerical 
calculations of the Eliashberg equations in the normal state [119] give very interesting 
behavior of the density of states N{(o), where a strong renormalization of N{(o) is 
present, but which is absent in the standard theory of the isotropic EPI. First, N{(0 = 
0) > Ni)are{(0 = 0), where Ni,are{(i> = 0) is the density of states in the absence of the EPI 
- see Fig. 21. 

There is a "pseudogap"-like feature in the region (^2/5) < (o <Q. where N{(o) < 
Nhareifo). The "pseudogap" feature disappears at T comparable with the phonon energy 
Q.. Note, that the usual isotropic EPI does not renormalize the density of states in the 
normal state, i.e. N{o)) = Nbareifo)- As a consequence of the pseudogap behavior of 
N{0)) the transport properties are very peculiar. For instance, the resistivity p{T) is 
linear in T starting at very low temperatures, i.e. p{T) ~ T for {D./30) < T and extends 
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FIGURE 20. The critical temperature Tc [K] of d-wave superconductor as a function of the experimental 
parameter p,/ ac[K], p, is the residual resistivity and a is defined in the text. The case /3 = corresponds 
to the prediction of the standard d-wave theory with isotropic scattering [140]. The experimental data 
[141] are given by crosses - YBa2{Cui-xZnx)-iO-]_s, and circles - Yi-yVvyBa2Cu2,Oqs - from [142]. 



up to several - as it is seen in Fig. 22. The dynamical conductivity (Ti(ft)) shows the 
(extended) Drude-like behavior with the Drude width F^^ ~ T, for (O <T - see Fig. 22. 
The above numbered properties are in a qualitative agreement with experimental results 
in HTSC oxides, as discussed in Section 2. 

In this extreme forward scattering peak limit one can calculate T^. In leading order 
w.r.t. (Q./Tc « 1) one has 

Tm^N{0)Vep = XN{0)/4, (97) 

where A = N{0)Vep. In that case the maximal superconducting gap is given by Aq = 2Tc 
which is reached on the Fermi surface, while away from it the gap decreases, i.e. 

Afc = Ao^l-(<^;t/Ao)2. (98) 

The expression for T^ tells us that it can be large even for A < 0. 1, since in HTSC oxides 
the bare density of states is Nbarei^) ~ Istates/eV. It is apparent that in this order there 
is no isotope effect, i.e. a = 0. We stress that such an extreme limit is never realized in 
nature, but for the self-energy it is a good starting point, since the effects of the finite 
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FIGURE 21. The density of states N{(o) in the FSP model for the EPI. with the dimensionless coupling 
/(= Vep/kQ) = 0.1 for various?(= TtT/Q.). From [119]. 



width (kc) of a^F(k, k', ©), whenever kc<^kp, change mainly the quantitative picture - 
see [1 19]. In case when kcVp <^ ^ the reduction of is given by 



K{2>)kcVF , 



cO 



(99) 



Very interesting calculations in the more realistic FSP model with the finite width 

kc, but kc ^ kf, were done in [134], where the FSP theory for the EPI and the 
SFI theory (based on spin-fluctuation mechanism of pairing) were compared. For in- 
stance, the FSP theory can explain the appreciable increase of the anisotropy ratio 
R = A{7t, 0) /A(;r/2, 7t/2) when T Tc, while the SFI is unable. Furthermore, the FSP 
theory of the EPI can explain the pronounced orthorhombic (a ^ b) effect in YBCO on 
the gap ratio Aa/Aj,, penetration depth anisotropy ^a/^b supercurrent ratio in the 
c-axis Pb — YBCO junction. On the other hand, the SFI theory is ineffective, since it 
predicts at least one order of magnitude smaller effects - [134], [2]. 




FIGURE 22. p{T) (upper part) and C7i(aj) in the FSP model for the EPI. with the dimensionless 
coupling /(= Vep/tcCI) = 0.1 for variousr(= nT/Q.). From [119]. 

6.3. Nonadiabatic corrections of Tc 

HTSC oxides are characterized not only by strong correlations but also by relatively 
small Fermi energy Ef, which is not much larger than the characteristic (maximal) 
phonon frequency CO^^^, i.e. ~ 0.1 — 0.3 eV, (O^^^ ~ 80 meV. The situation is 
even more pronounced in fullerene compounds A3C60, with Tc = 20 — 35 K, where 



Ep — 0.2 eV and O)^^^ ~ 0.16 eV . This fact implies a possible breakdown of the 
Migdal's theorem [139], [9], which asserts that the relevant vertex corrections due to 
the £■ — P interaction are small if {(Oo/Ep) <^ 1. In that respect a comparison of the 
intercalated graphite KC^ and the fuUerene A3C60 compounds, given in [120], is very 
instructive, because both compounds have a number of similar properties. However, the 
main difference in these systems lies in the ratio (Od/Ep, since (cod/Ef) <^ 1 in KC^, 
while it is rather large {(Oo/Ep) ~ 1 in AsCgo- Due to the appreciable magnitude of 
(On/Ep in the fuUerene compounds and in HTSC oxides it is necessary to correct the 
Migdal-Eliashberg theory by vertex corrections due to the EPI. It is well-known that 
these vertex corrections lower in systems with isotropic EPI. However, the vertex 
corrections in systems with the forward scattering peak and with the cut-off Qc «kp 
the increases of Tc appreciable. The calculations by the Pietronero group [120] gave 
two important results: (1) there is a drastic increase of T^ by lowering Qc = qd^kp, for 
instance Tc{Qc = 0.1) ~ ^Tc{Qc = 1); (2) Even small values of A < 1 can give large T^. 
The latter results open a new possibility in reaching high T^ in systems with appreciable 
ratio (Od/Ep and with the forward scattering peak. The difference between the Migdal- 
Eliashberg and non-Migdal theories can be explained qualitatively in the framework of 
an approximative McMillan formula for T^ (for not too large A) which reads 

TcP^ {(o)e-^l^^'^*\ (100) 
The Migdal-Eliashberg theory predicts 

while the non-Migdal theory [120] gives 

A^A(1+A). (102) 

For instance Tc ^ 100 in HTSC oxides can be explained by the Migdal-Eliashberg 
theory for A ~ 2, while in the non-Migdal theory much smaller coupling constant is 
needed, i.e. A ~ 0.5 as it is seen in Fig. 23. The pioneering approach done in [120] 
deserves more attention in the future. 



6.4. Pseudogap behavior in the FSP model for the EPI 

In this review we did not discuss a number of interesting topics such as the possible 
existence of stripes, the properties of the pseudogap state, etc.. There is a believe that 
the understanding of these properties might give some hints for pairing mechanism in 
HTSC oxides. Especially, the pseudogap (PG) problem is a very intriguing one and is not 
surprising at all, that a number of theoretical approaches were proposed for explaining 
the PG. We are not going to discuss it but only quote some of them. first one is based 
on the assumption that the PG phase represents pre-formed pairs [143], and the true 
critical temperature T^ is smaller than the mean-field one T^^. In the region Tt<T<T^^ 
pre-formed pairs exist giving rise to the dip in the density of states N{(0). This approach 
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FIGURE 23. The approximative analytic (solid lines) and numerical (circles) solution of Tc{X) in the 
first nonadiabatic approximation for various cutoff Qc and for m = {(Od/Ef) = 0.2. From [120]. 



is physically plausible having in mind that HTSC oxides are characterized by the short 
coherence length and quasi-two dimensionality. From the experimental side there are 
some supports. For instance, the specific heat measurements [144] point to the non-mean 
field character of the superconducting phase transition, particularly for the underdoped 
systems. As we have already mentioned in the Introduction the ARPES measurements 
show that the PG has a d-wave like form 

Apg(k) w IS.pg^Q{co%kx — cosfcy) (103) 

(like the superconducting gap) and Apg o increases by lowering doping. The second 
approach assumes that the PG is due to a competing order, but usually without the 
long-range order, such as due to "spin-density wave" alias for strong antiferromagnetic 
fluctuations [145], [146]. There are other approaches which are based on the RVB and 
orbital current model, d-wave order, etc., but we shall not discussed it here. 

However, the FSP theory, which predicts the long-range force due to the renormaliza- 
tion of the EPI by strong correlations, opens an additional possibility for the PG. As we 



discussed in Section 6.2, due to the forward scattering peak the critical temperature has 
a non-BCS dependence, i.e. T^^ — Vep/^- However, this is the mean-field value, which 
is inevitably reduced by phase and internal Cooper pair fluctuations present in systems 
with long-range attractive forces, i.e. with the forward scattering peak. 

The interesting problem of fluctuations in systems with long-range attractive forces 
was recently studied in [147]. It was shown there, that such a long -ranged superconduc- 
tor exhibits a class of fluctuations in which the internal structure of the Cooper pair is 
soft. This leads to a "pseudogap" behavior in which the actual transition temperature 
Tf is greatly depressed from its mean-field value . We stress that these fluctuations 
are not the standard phase fluctuations in superconductors. Since the problem is very 
interesting and deserve much more attention in the future we discuss it here briefly. In 
the following the weak coupling limit is assumed, where the pairing Hamiltonian has 
the form 

// = J: /rfxVAt(x)^(p)V/^(x) 
a •' 

-J JxifxV(x-x')VA|(x)v^|(x')VA|(x')v/|(x). (104) 
In the MFA the order parameter A(x, x') is given by 

A(x,xO = y(x-xO(v^i(xOv/tW)- (105) 

A(x, x') depends in fact on the internal coordinate r = x — x' and the center of mass 
R = (x + x')/2, i.e. A(x,x') = A(r, R). In usual superconductors with the short-range 
pairing potential ys^(x — x') pa Vo5{x — x') one has A(r,R) = A(R) and therefore there 
are practically the spatial (R) fluctuations of the order parameter, only. In the case of 
long-range pairing potential there are additional fluctuations of the internal (r) degrees 
of freedom. In the following we sketch the analysis given in [147]. 

When the range of the pairing potential is large, i.e. Tc > t, (the superconducting 
coherence length), fluctuations of the internal Cooper wave-function are important since 
they give rise to a tremendous reduction of the mean-field quantities. In order to make 
the physics of internal wave-function fluctuations we study much simpler Hamiltonian 
the so called reduced BCS Hamiltonian, 

H = Y, ^^Aa^^o - Vk_k/4^cl^^c_k/|Ck/^ . (106) 
kc7 k,k' 

Since we shall study excitations around the ground state we assume that there are no 
unpaired electrons which allows us to study the problem in the pseudo-spin Hamiltonian 
[148] 

k<j ^ k,k' 

= E2^k5L - E y^MSlSi,+Sl,Si), (107) 
kc7 k,k' 



where the pseudo-spin 1/2 operators 5^^, S^^ = {S^-^y are given by 

We see that Eq.(107) belongs to the class of the Heisenberg ferromagnetic (Vjj-k' > 0) 
Hamiltonian formulated on the lattice in the Brillouin zone. The mean-field approxima- 
tion (MFA) for this Hamiltonian is given by 

//mfa = -l^hkSk (109) 

k 

with the mean-field hk given by 

hk = -2<^kZ+J^Vk-k'(5^'X+5];,y), (110) 

k' 

where x,y and z are unit vectors. Since x- and y-axis are equivalent one can searches 
hk in the form hk = — 2^k2;+2AkX, where the order parameter Ak is the solution of the 
equation 

Ak = i:Vk_k'(5^k') = I^k-k'^^^tanh^, (111) 
k' k' k 

with Ek = + 

In the case of short-range BCS-like forces Vbcs{^ — x') ~ VqS (x — x') one has Vk-k' — 
Vo for all momenta. This "long-range force" in the momentum space it is the "long-range 
force" means that the MFA is good approximation with the standard BCS solution of 
Eq.(lll). 

For the long-range attractive forces the function Vk-k' is peaked at | k — k' | = 0, 
for instance in the extreme forward scattering peak case (see Section 6.2) one has 
Vk-k' = ^o^(k — k')- the following we analyze s-wave pairing only where the solution 
of Eq.(lll) gives T,^^ = Vq/A and Aq = IT^^ . (Note, that in the BCS case one has 
Aq = \.16T^^ .). The coherence length is defined by ^ = vp/n/^. The important fact 
is that in the case of long-ranged superconductors the Heisenberg like Hamiltonian in 
the momentum space is short-ranged giving rise to low-laying spin-wave spectrum. 
The latter spectrum are in fact the low-energy bound states (excitons) which loosely 
correspond to the low-energy collective modes (in the true many-body theory based 
on Eq.(104)). This problem is studied in [147] for the long-range (but finite) potential 
V{r) = Vbexp{— r^/2r^} (its Fourier transform is Vk — (2;rr^)Vbexp{— A:^r^/2}) where 
it was found a large number Ncm ~ '^I^fT'cI (for r^O^ ) of the excitonic like collective 
modes at zero momentum. These excitonic modes lie between the ground state and 
the two particle continuum for © > 2Ao. Note, that since we assume that Aq <C Ep the 
system is far from the Bose-Einstein condensation limit. 

The above analysis is useful for physical understanding, but the fully many-body fluc- 
tuation problem, which is based on the Hamiltonian in Eq.(104), is studied in [147] 



where the Ginzburg-Landau (G-L) equation is derived for the long-ranged superconduc- 
tor. Due to the fluctuations of the internal wave-function the G-L free-energy functional 
F{A(R,k)} for the order parameter A (R,k) = / i/rA(R-r/2,R-Fr/2)exp{ /kr} has 
much more complicated form 

F{A(R,k)} = £ / dR{Ay, I A(R,k) \^ +5k | A(R,k) \^ 

k 

1 ,2 , 1 ,:...„.^,2^ 



where M = riVo and 



+— I dkA{R,k) \' +^ I aRA(R,k) 1^}, (112) 



^^^J_ tanh(j8W2) 



Vb 2^k 

1 j82v2sinh(j8^k/2) 



5k = ^r^f^-TTT^ (113) 



2mk 32^kCOsh^(j8^k/2) 

tan h(j8^k/2) ^ 

8^3 16^2,o3h3(j8^k/2) 

The term due to the partial derivative (9k is a direct consequence of the long-ranged 
pairing potential, and it describes of fluctuations of the internal Cooper wave-function. 
The effect of these fluctuations, described by the free-energy functional in Eq.(112), is 
studied in the Hartree-Fock approximation in the limit ^ , where it is found the 
large reduction of the mean-field critical temperature 

Tc-^r^y (114) 

The latter result means that Tc in the long-ranged superconductors is controlled by 
thermal fluctuations of collective modes which is in contrast with the short-range (BCS- 
like) superconductivity. In the temperature interval Tc < T < T^^ the system is in 
the pseudogap regime where the electrons are paired but there is no long-range phase 
coherence. The latter sets in only at T < Tc. We shall not further discuss this interesting 
approach but only stress that it can be generalized by including the repulsive interaction 
due to spin fluctuations, what shall be discussed elsewhere. 

In conclusion, the forward scattering peak in the EPI gives rise to the long-ranged su- 
perconductivity in which the soft excitonic modes of the internal Cooper wave function 
reduce Tc strongly. In the region Tc<T < T^^ the pseudogap (PG) phase is realized. In 
this approach the PG has the same symmetry as the superconducting gap. 



7. ELECTRON-PHONON INTERACTION VS 
SPIN-FLUCTUATIONS 



7.1. Interaction via spin fluctuations (SFI) and pairing 

At present one of the possible candidates in explaining experimental results in HTSC 
oxides appears to be the theory based on the spin fluctuation pairing mechanism - the 
SFI theory. The latter is usually described by the single band Hubbard model, or on the 
phenomenological level by the postulated form of the self-energy (written below)[29], 
[136], [135], [83], [84], [137]. In the approach of Pines- school to the SFI the effective 
potential Ve//(k, to) (see Eq.(5) in Sections 2.) depends on the imaginary part of the 
spin susceptibility Im;^(k — k', co) (co real). According to this school, the shape and the 
magnitude of Im;^(q, (o), which is peaked at Q = {7t,n), plays an important role in 
obtaining in this mechanism. There are two phenomenological approaches, which can 
be theoretically justified in a very weak coupling limit gsf <C 1 only, where Im;^(q, co) 
is inferred from different experiments: 

(1) From NMR experiments at very low (D - the MMP model [29], [136], [135], where 
IfnXuMP is modelled by 



7.,„„(,,<.+,-o-) = ^ „^^.|^_^^°^(,^,^^).,. e«-- 1 « I), (115) 

with the frequency cutoff (O^^^ = 400 meV . They fit the NMR experiments by assuming 
very large value for Xq ~ (30 — AO)xq ~ 100 eV^^ . From Fig. 24 it is seen that the 
imaginary susceptibility is peaked at low frequency (Opeak ~ 5 — 10 meV . 

(2) From the neutron scattering experiments [83], [84]) - the RULN model, where 
I^Xruln is modelled by 

, 1 , 
ImXRULNiq, (O + j0+) = C[—— r]^ x 

'1-|-/0[COS^^-|-COS^3,] 

X 3(r + 5)ft) ©fffl^t^^- I ffl I) ril6^ 

1.50)2-60 I CO I +900 + 3(r + 5)2^^'^- I 1^' ^^^^^ 

where (O^^^^ = 100 meV, Jq = 0.3, C = 0.19 eV^^ with T and (O measured in meV. 
From Fig. 25 it is seen that ImXRULN is peaked around 30 meV, which is much larger 
than in the MMP model. 

By knowing Imx one can calculate the effective pairing potential Ve//(k, O)) from 
Eq.(5) and the spectral function for the d-wave pairing ocJF(o) 

Here, }d(k) = cos^;;^ — cos^y is the d — wave pairing function (A(k, (o) ^ A{(o)Yj{k)). 
The bracket means an average over the Fermi surface. The spectral function a^F{(o) for 
two models is shown in Fig. 26, where it is seen that ajF((o)^^^^ is much narrower 




FIGURE 24. Spectral function lmx{, (o) for the MMP model of spin-mediated interactions \mx{, (o) 
in YBa2CusO-j_s. The spectral function is calculated at = {7t,7i) and T = K (solid Une), 100 K 
(long-dashed line), 200 K (short dashed line), and 300 K (dot-dashed line). Inset: experimental data of 
YBaiCusOfy.e atT = 100K- the hne is to guide the eye. From [83], [84]. 



function than ajF(ft))*'^^. The latter is peaked almost at the same (O as ajF(o))^'^^", 
while ajF{co)^^^ is much broader than ajF{a)Y^^^ . 

Due to different shapes of the susceptibility and of a^F{(o) in these two approaches 
the calculated (from Eliashberg equations) critical temperatures are also very different. 
Since the MMP spectral function is much broader than the RULN one it turns out that 
rj,{MMP) j-each 100 ^ for rather large value of gsF ^ 0.64 eV , while yj^^^^) saturates 
already at 50 K even for gsF ^ 1 • From the physical pint of view the RULN model is 
more plausible than the MMP one, since the former is based on the neutron scattering 
measurements which comprise much larger frequencies than the l^^MR measurements. 
Note, that a valid model for HTSC oxides must be able to explain the high values of Tc 
(which needs Ksf{= 2 / (ajF((y) / (o)da)) ~ 2) and the resistivity p [T] (and its slope p' 

with small Xtr ~ 0.6). It turns out that jj^^^^^ can fit Tc ~ 100 K on the expense of large 
coupling g[^^^^ ~ 0.64 eV and X^^^^^ ~ 2.5. However, the value g[^^^'' ~ 0.64 eV 
gives much larger value for p(r) and p' than the experiments do. On the other hand if 
one fits p{T) and p' with the MMP model one gets very small Tc <1 K, thus making 
the MMP model ineffective in HTSC oxides. 

In the physically more plausible RULN model jj:^^^^^ saturates at 50 K even for 

sff"^^^ ^ 1 • chooses an appropriate value for g^^^'^ to fit p(T) and p' one 
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(0 (meV) 

FIGURE 25. Spectral function lmx{,(o) for the RULN model of spin-mediated interactions x(,fi)) 
in YBa2CuiO-i_s- The spectral function is calculated at = {7l,n) and T = K (solid line), 100 K 
(long-dashed line), 200 K (short dashed line), and 300 K (dot-dashed line). Inset: experimental data of 
YBa2Cu306.6 at r = 100 /sT - the line is to guide the eye. From [83], [84]. 
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FIGURE 26. The spectral function a^F{(o) in the d — wave pairing channel for the MMP and RU LN 
model at different temperature - from [137]. 



gets j^^^^^^ ^ 1 K. This analysis gives a convincing evidence that the existing SFI 
theories are ineffective in HTSC oxides. 

We stress again, that the large effective coupling constant, assumed in the SFI theories, 
gsf ~ 0.64 eV , is difficult to justify theoretically (if at all). By analyzing theoretically 

the possible strength of the coupling constant g^f, in both the weak (N{0)U ^1) and 
strong {N{0)U ^1) coupling limit, one obtains that g^f < 0.2 eV and X^f < 0.2 (note 

Xsf ~ g^jr), which means that Xsf ^ xj,'^'^^^ and rj'^-^"* <^ ji'^^P) analysis is 

supported by the recent theoretical results where gsf < 0.2 eV is extracted from the 
calculation : 1. of the width of the magnetic resonance peak at 41 meV [68]; 2. of the 
small magnetic moment (jl < 0. 1 Hb) in the antiferromagnetic order, which coexists 
with superconductivity in La2-xSrxCuO^ [69]. 



7.2. Are the EPI and SFI compatible in d-wave pairing? 

The phenomenological SFI theories (MMP, RULN, FLEX approximation [149]] be- 
came popular because they can produce d-wave pairing, due to the repulsive character of 
spin-fluctuations (in the momentum space.) which are peaked in the backward {q ~ Q) 
scattering. However, as we have argued in previous Sections, a number of experiments 
point to a large EPI with X > 1 . On the other hand if one assumes that the EPI is mo- 
mentum independent (isotropic), like in the standard Migdal-Eliashberg theory, then it 
is strongly pair-breaking for d-wave pairing. So, if one assumes (for the moment) that 

is f) 

superconductivity in HTSC oxides is due to the SFI with Xsf ~ 2, then in that case Tc 
would be drastically reduced (to almost zero) by the isotropic EPI even for moderate 
Xep ~ 1. The latter was shown in [138] where the Eliashberg equations are solved for 
the SFI treated in the FLEX approximation [149] and the EPI in the Einstein model 
with various momentum dependent VEp(k, (d). This result means, that if the SFI would 
be the basic pairing interaction in the presence of the isotropic and momentum inde- 
pendent EPI, then in order to reach ~ 100 ^ the bare critical temperature should be 

(s f) 

Tc ~ (600 — 700) K, which needs unrealistically large Xsf. This is not only highly 
improbable, but would give enormous large resistivity and its slope, in contrast to exper- 
iments. The similar results were obtained in [1 1]. 

The calculations in [138] show that the SFI interaction is dominant in (d-wave) pairing 
if some strong constraints are realized, such as: (1) very large SFI coupling constant 
Xsf ^2; (2) a strong forward scattering peak in the EPI with small EPI coupling A <C 1. 
Both these conditions are incompatible with experiments and theoretical analysis - see 
also Section 2.1. The way out from this controversy is that the EPI with the forward 
scattering peak is inevitably dominant interaction in the quasiparticle scattering and 
pairing in HTSC oxides. As we already discussed in Section 5. the forward scattering 
peak in the EPI gives rise to the large coupling constant in the d-wave channel, which is 
of the order of the one in the s-wave channel in the range of doping around the optimal 
one, i.e. Xsp^d ~ ^ep,s- This means that the residual Coulomb repulsion (by including 
also the SFI with the backward scattering peak (BSP)) with Xc < Xgp^d triggers d-wave 
pairing. 



8. IS THERE HIGH-TEMPERATURE SUPERCONDUCTIVITY IN 
THE HUBBARD AND T-J MODEL? 



8.1. Hubbard model 

There are a number of papers dealing with numerical calculations, such as Monte 
Carlo, exact Lanczos diagonalization, in the 2-dimensional (2D) single-band and three- 
band Hubbard model. One can say that the single-band Hubbard model does very well 
in describing the magnetic properties of HTSC oxides. Concerning the existence of su- 
perconductivity the situation is not definitely resolved. So far the calculations are done 
on finite clusters and rather high temperatures T > O.U [150], [7] which show no ten- 
dency to superconductivity. It is worth of mentioning that most of these calculations 
deal with the pairing susceptibility - see Eq.(126) below, defined in terms of the bare 
electron operators c^athbtka in Eq.(128). Since superconducting pairing is realized on 
quasiparticles with the weight z < 1 there is a last hope that the accuracy of the present 
numerical calculations is not sufficient to pick up the suppressed pairing susceptibility. 
In that respect a very important approach to the problem of superconductivity (in any 
microscopic model), which is formulated without using any order parameter, was given 
by the Scalapino group [151]. The method of calculations is based on the two most im- 
portant hallmarks of superconductivity: (i) ideal diamagnetism (the Meissner effect) and 
(ii) ideal conductivity. In that respect they study the superfluid density Dg (proportional 
to A^^, A is the penetration depth) and the Drude weight D in the single-band n.n. (near- 
est neighbors) Hubbard model. The dynamical conductivity along the x-axis is given 

I (O + IO 

Here, (— 7i) = {—T) /2 where T is the kinetic energy in the n.n. tight-binding model - 
see Eq.(18), where the current-current response function Axx{q, CO) is obtained from 



AU^,i(Om) = ^ [^dze'"^'{jP{q,z)jP{-q,0)), (119) 
i\ Jo 



with (Om — iKmT, by the standard analytic continuation ico^ (0 + 15 and 



7^ (q, ^) = ^■^E^"'''VLcT^/a - clci+,a). (120) 
l,a 

In the pure Hubbard model Oxx{(o) contains the delta function contribution 

O^^{C0)=Dd{C0) + Oreg{C0), (121) 

where the Drude weight D{= (n/m)*, which measures the ratio of the density of the 
mobile charge carriers to their mass, is defined by 

^ = i-T,) - A«(q = 0, ft) ^ 0). (122) 



The Meissner effect is the current response to a static and transverse gauge potential 
q • A(q, CO — 0) — 0. In the small q limit one has 

(7a(q^ 0)) = -^{5ap -qaqp/q^)Ap (123) 

where Z)^(= (ns/m)* 

Ds = i-T^) - h^{qx = 0,^3, ^ 0, © = 0). (124) 

Based on the above definitions of D and we can study various phases of the system: 
(1) D=D^=0 - an isolator, (2) Dt^O and Dj=0 - a nonsuperconducting metal; (3) 7^0, 
D^O - a superconducting metal. 

The numerical Monte Carlo calculation in the repulsive Hubbard model (U=4t>0) 
[151] on an 8x8 lattice show that Dj=0 and Dt^O in a broad range of the filling 
0.5<n<0.9 and for T > O.lt. This means that there is no tendency to high-temperature 
superconductivity in the single-band Hubbard model. This conclusion is supported by 
the projector-QMC calculations [152] for the quarter filling case n=0.5 and at T=0. 
That these results (Ds = 0) are not a finite size effect confirm the calculations on the 
attractive Hubbard model (U=-4t<0), also on an 8 x 8 lattice, where the clear tendency to 
superconductivity is found, since Dy 7^0, D^^O already at T<0.2t. 

The paper [151] is of great importance numerical studies of pairing in model systems, 
not only because it hints on the absence of superconductivity in the repulsive Hubbard 
model, but also because of the following two reasons: (1) It uses the general and unbiased 
criterion for superconductivity, which is independent on the type of the pairing ampli- 
tude; (2) It shows that the attractive interaction is more favorable for (high-temperature) 
superconductivity than the repulsion. 



8.2. t-J model 

The SFI phenomenological approaches root on their basic t — / Hamiltonian Eq.(61). 

On can put a legitimate question - is there superconductivity in the t-J model? In the 
past there were various approaches confronting with this important problem. In spite of 
a number of controversial statements it seems that the results converge to the unique 
answer - there is no superconductivity with appreciable Tc. If superconductivity exists 
Tc is very low. As the strong support of this claim serve the recent calculations based on 
the high-temperature expansion in the t-J-V model [70], 

^=-i:^0-^f°^f + E m-Sj) + {V-^-)mnj)]. (125) 
i,j,a i,j=n.n. 

Here, the V-term mimic the screened Coulomb interaction which is always present in 
metals, where one expects that V > J.Ia [70] it was calculated the uniform susceptibility 
for the superconducting pairing 

Xsc =ljj^ dTiiW^Osce-'^'olc)) (126) 
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FIGURE 27. Superconducting (d-wave) susceptibility %{d — Sc) (T) for t ~ 2J. Pairing correlations are 
already weak for V = 0, dashes, and decreases further by increasing V-n.n. repulsion (dot-dash). From 
[70] 



05C=^£(Ar,r+x±Ar,r+y) (127) 
^ r 

^'j = C'T^ji+^7T^U- (128) 
where + sign holds for the s-Sc and the - sign for d-Sc. 

In the physical region of parameters t > J both Xs-sc and Xd-sc are small and further 
decrease by decreasing T. For rather small V = J/4, which is even much smaller than 
expected, the superconducting susceptibilities are drastically decreased as it is seen in 
Fig. 27. This means that in the more realistic models for HTSC oxides, such as the t-J-V, 
there is no tendency to high-temperature superconductivity. If superconductivity exists 
at all its Tc must be very low. 



9. SUMMARY AND CONCLUSIONS 

A number of experiments, such as optics {IR and Raman), transport, tunnelling, ARPES, 
neutron scattering, give convincing evidence that the electron-phonon interaction (EPI) 
in HTSC oxides is sufficiently strong and contributes to pairing. These experiments give 
also evidence for the presence of strong correlations which modify the EPI not only 
quantitatively but also qualitatively. The most spectacular result of the EPI theory in 
strongly correlated systems is the appearance of \hQ forward scattering peak in the EPI, 



as well as in other charge scattering processes such as the residual Coulomb interaction 
and scattering on non-magnetic impurities - the FSP theory [17], [18], [19], [20] 
[2]. The forward scattering peak is especially pronounced at lower doping 5. This 
fundamental result allows us to resolve a number of experimental facts which can not 
be explained by the old theory based on th isotropic Migdal-Eliashberg equations for 
the EPI. The most important predictions of the FSP theory of the EPI and other charge 
scattering processes are: (1) the transport coupling constant Xtr (entering the resistivity, 
p ~ XtrT) is much smaller than the pairing one A, i.e. Xt,- < A/3; (2) the strength of 
pairing in HTSC oxides is basically due to the EPI, while the residual Coulomb repulsion 
(including spin-fluctuations) triggers d-wave pairing; (3) d-wave pairing is very robust 
in the presence of non-magnetic impurities; (4) the nodal kink in the quasiparticle 
spectrum is unshifted in the superconducting state, while the anti-nodal singularity is 
shifted. 

We stress the following two facts coming from the theoretical analysis: (i) the forward 
scattering peak in the EPI of strongly correlated systems is a general phenomenon by 
affecting electronic coupling to all phonons; (ii) the existence of the forward scattering 
peak in the EPI is confirmed numerically by the Monte Carlo calculations for the 
Hubbard-Holstein model with finite U [21], by exact diagonalization [124], as well as 
by some other methods [22]. 

Tunnelling experiments and ARPES measurements of the real part of the self-energy 
give evidence that the EPI coupling constant A > 1. At present there are no reliable 
microscopic calculations of A in HTSC oxides, which properly include (a) the ionic- 
metallic coupling (due to the long-range Madelung energy) and covalent coupling and 
(b) strong electronic correlations. 

In the last several years a large number of published papers were devoted to the 
study of spin-fluctuation (SFI) interaction as a mechanism of pairing in HTSC oxides. 
In spite of many efforts and well financed theoretical projects (headed by the greatest 
authorities in the field), which have opened some new research directions in the theory 
of electron magnetism, there is no theoretical evidence for the effectiveness of the non- 
phononic mechanism of pairing. Until now superconductivity could not been proved in 
the repulsive single-band Hubbard model as well as in its derivative the t — J model. 
Just opposite, quite recent numerical calculations [70] show in a convincing way, that 
there is no high-temperature superconductivity in the t-J model. If it exists its T^ is 
extremely low. Finally, the numerical calculations in the Hubbard model [151] show that 
the repulsive Hubbard interaction is unfavorable for high-temperature superconductivity, 
contrary to the attractive interaction which favors it. 

The explanation of the high critical temperature in HTSC oxides should be searched in 
the electron-phonon interaction which is renormalized by strong electronic correlations. 

To conclude, one can not avoid unavoidable. 
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10. APPENDIX: DERIVATION OF THE T-J MODEL 

10.1. Hubbard model for finite U in terms of Hubbard operators 

For simplicity we study the nearest neighbor (n.n.) Hubbard model [153] 

H=-t 4a^n(T + f/J^«mT«m^ (129) 

m^na m 

where the operator c\i(j creates an electron at the m-th site with the spin projection o. 

The Hilbert space at the given lattice site contains four states { | a >^^\ >, | 2 > 
,|t>7|i>}- Let us introduce the Hubbard projection operators Z"^ ; oc, j8 = 0,2,(7 
(where a =t (+),a=i (-)) 

=1 a >< j8 I . (130) 

They fulfill the projection properties 

= (131) 
and rather "ugly" (anti)commutation algebra 

Xfxf±xfxf = 5,j{5^yXf ± dsaX^^. (132) 
The completeness relation of the Hilbert space reads 

The Hubbard operators describe the composite object. There is a connection between 
QaandX«^ (if(7=T==^a=i) 

= Xf" + aZp; 4 = + oXf^ (134) 



n/=l-XPO+zP (135) 
g+ — fit fi., _ v+- — rc-^t _ rv-+^t 



. =£lfii=xr = (S7y = (xry 



and vice versa 



X''^ = cl{l-ns) (137) 
Z^*^ = n^(l-/ia); X'^* = 4ca (138) 

X^=(l-n^){l-ni) (139) 

Z2^ = (7ctn,,; Z20 = (7cjcc (140) 
X^^ = n^ni (141) 
The Hubbard Hamiltonian H — Hi+ H12 + //q in terms of X"'^ is given by 

Hi = -tY, {X^'^X^'' + (142) 
= 1^ (7(X;'^°Zf 2 +xpZj"^) (143) 

Ho = uY,Xl^ (144) 

The first term in Hi describes the motion of single electron in the lower (L) Hubbard 
band, while the second term describes the motion of the doubly occupied electrons from 
j-th to the i-th side in the upper (U) Hubbard band. The term Hq is the repulsive energy 
of two electrons on the i-th site. Finally, H12 connects the two (lower and upper) bands. 



10.2. Effective Hamiltonian f or U»t 

There are various ways to obtain the effective Hamiltonian in the case U » t. 
Because of its generality and simplicity we use here the canonical transformation method 
[153], where the operator S mixes lower and upper band. Under the action of S the 
Hamiltonian is transformed into Heff = e^He~^ 

Heff = H+[S,H] + ^[S, [S,H]] + .. (145) 

with S in the form 



Now, we choose K" so that all first-order in t processes between the L- and U-band 
disappear from Hgff, i.e. one has 

Hn+[S.H2]=Q. (147) 
The solution of Eq.(147) isK = —t/U, and Heff reads 

Heff = -t Y^X-^^X^^ + H2,s 
ija 

+jY^{SiSj-hinj)+H2, (148) 
ija 

where / = 2t^/U is the exchange energy. 

The term H2 describes motion of "doublons" in the U-band 

H2 = f/£xp -tJ^X^^Xp, (149) 
i ija 

while H^s describes the three-sites hopping. 

H3s = ^ (Xf^xrxf^ -Xf^X^X^^). (150) 
ijla 

Usually this term is neglected in the t-J model. 

However, it may have tremendous effect on superconductivity by strongly suppressing 
it [154]. By projecting Hgff onto the lower Hubbard band one gets the famous t-J model 
Hamiltonian Hfj = PHg/fP 

H,j = -t J^Zf + +/ J^(S,-Sy - hafij) 

ija ija 

= £xf + i- Y^iXf^^xf -X^^xf). (151) 

ija ija 

Before we are going to discuss some representations for non-canonical operators X^^ in 
terms of bosons and fermions let us stress that the so called "spin" operators 5*, 5^ do not 
describe correctly the electron spin. Although they satisfy the correct spin-commutation 
relations 

[S],Sj]^±5ijSf, (152) 
they describe a particle with spins 5 = 0, 1/2 at the lattice site, since S? fulfills 

^] = \ni^\. (153) 



Since X"^ obey the non-canonical ( "ugly") algebra the question is how to treat the 
Hamiltonian Htj ? In Section 4.-5. we have shown that one can study directly with 
these operators by using the functional technique and 1/N expansion for the self-energy. 
However, there are very popular approaches which represent Z"^ in terms of bosons 
and fermions with canonical commutation relations. 

Slave boson (SB) method. Here, one introduces the fermion with spin (spinon) i^^ 
and the boson without spin (holon) 5,- operators, where X^^ = F(jB'^. The constraint 
on the SB Hilbert space (completeness), at the given lattice site, is given by B^B + 
Lc7 ^aPcF = 1 and the t- J Hamiltonian reads 

HtJ = -tY^FicFjcBiB] + ^ £ FlFjaFl,Fia. (154) 
ija ijoa' 

We stress that the constraint strongly limits the SB Hilbert space of bosons and fermions 
(at the given lattice site) which effectively means their strong interaction. In that respect 
any uncontrollable decoupling in the SB method (as in some RVB approaches) leads to 
a spin-charge (spinon-holon) separation, which is not realized in 2D and 3D systems. In 
order to correct this one introduces the so called gauge fields which keep the spin and 
charge together. We already discussed the difficulties of the SB method in studying the 
electron-phonon interaction. 

Slave fermion (SF) method. In the SF method the boson has spin and fermion not, 
i.e. X^^ = b]jF with the constraint on the Hilbert space F^F + Ler^CT^cy = 1- 

Spin fermion method. Here, the real fermion with the "spin" S is represented via 
the auxiliary fermion F^ and spin s by c|c| -|- cjcj = 1 — F^F and S = s(l — F^F). The 
t-J Hamiltonian is rather complicated 

Hu = 2tY,F,^Fj{siSj + ^) 

ij 

+jY^{l-F^^Fi){siSj-\){l-FjFj). (155) 

ij 

Usually this method is used for analyzing motion of single hole in the half-filled system 
where the antiferromagnetic order is realized. 
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